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LIGHT  SCATTERING  FROM  BUBBLES  IN  LIQUIDS 
Abstract 

by  Dean  Scott  Langley,  Ph.D. 
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Chair:  Philip  L.  Marston 

The  scattering  of  laser  light  from  bubbles  In  liquids  was 
photographed,  and  angular  features  in  the  far-zone  intensity  are 
compared  with  models.  Diffraction  effects  are  prominent  in  three 
regions  of  the  scattering:  near  the  critical  angle,  and  in  the 
forward  and  backward  directions.  The  forward  and  backward  scattered 
light  exhibits  the  intensity  enhancement  known  as  the  glory.  Simple 
geometrical  optics  is  unable  to  approximate  the  scattering  pattern  in 
the  diffraction  regions.  Physical-optics  models  are  presented  for 
each  case  and  shown  to  be  in  agreement  with  the  observations.  These 
models  give  insight  into  the  optical  properties  of  scatterers  whose 
refractive  index  is  less  than  that  of  their  surroundings. 
Comparisons  also  include  the  results  of  Mie  theory  computations  for 
the  critical-angle  and  forward  regions.  Some  applications  of  the 
results  to  the  sizing  and  detection  of  bubbles  are  discussed. 
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CHAPTER  1 


INTRODUCTION 


1 .1  Overview 

This  dissertation  is  concerned  with  the  light  scattering 
properties  of  bubbles.  A  bubble  is  distinguished  from  a  droplilce 
scatterer  by  a  refractive  index  n^  which  is  smaller  than  that  of 
the  surrounding  medium  nQ;  thus,  the  relative  index  m  »  n^/n^  <  ^ 
for  bubbles  while  ra  >  1  for  drops.  The  literature  on  light 
scattering  contains  comparatively  little  Information  on  the 
scattering  properties  of  bubbles,  but  significant  differences  from 
the  scattering  by  drops  can  be  observed.  The  purpose  of  this  thesis 
is  to  present  measurements  and  models  of  the  scattering  from  bubbles, 
emphasizing  certain  angular-  regions  which  cannot  be  adequately 
described  by  simple  geometrical  optics.  There  are  three  such 
regions,  and  these  constitute  the  major  subdivisions  of  che  present 
research:  the  forward  and  baclcward  directions,  where  exceptionally 
strong  (glory)  scattering  is  present,  and  at  the  critical  angle  where 
the  transition  from  partial  to  total  reflection  of  rays  at  che 
bubble's  surface  occurs.  In  each  of  these  regions,  a  correct 
description  of  the  scattering  requires  that  diffraction  effects  be 
included,  even  for  bubbles  much  larger  than  the  wavelength  of  light. 
Features  in  the  scattered  light  can  provide  useful  information  about 
the  bubble,  while  a  successful  model  can  give  insight  into  the 
scattering  process. 
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The  problem  of  determining  the  pettern  of  scattered  light 
from  a  particle  of  known  shape,  size,  and  composition  is  generally  a 
difficult  one.  In  some  cases,  complete  solutions  of  Maxwell's 
equations  can  be  obtained,  giving  the  exact  amplitude  and  phase  of 
the  scattered  electromagnetic  field  everywhere.  Though  such 
solutions  are  certainly  valuable,  they  are  not  always  the  most 
desirable.  They  often  present  significant  computational 

difficulties,  and  the  effects  of  changes  in  characteristics  of  the 
particle  may  be  hard  to  assess.  In  many  cases,  approximate  solutions 
can  be  developed  having  a  limited  domain  of  applicability  but 
clarifying  the  important  physics  of  problems  within  that  domain.  Of 
course,  there  is  much  more  physical  insight  to  be  gained  when  an 
approximate  theory  can  be  compared  on  common  ground  with  an  exact 
one,  and  when  both  can  be  checked  against  experiment.  This  is  the 
approach  followed  here  when  possible.  Comparisons  of  models  and 
measurements  are  based  mainly  on  the  angular  locations  of  features  in 
the  scattered  intensity.  Throughout  this  thesis,  intensity  is  used 
to  denote  the  modulus  of  the  Poyntlng  vector. 

With  increasing  size,  bubbles  in  liquids  pass  from  spherical 

shapes  into  an  ellipsoidal  regime,  and  very  large  bubbles  take  on  a 

spherical-cap  form.^  The  bubbles  studied  in  this  work  were  much 

larger  than  the  wavelength  of  light,  but  small  enough  to  retain  the 

near-spherical  form.  The  lower  size  limit  removes  the  scattering 

2 

from  the  realm  of  the  Rayleigh  small- par  tide  approximation.  Nearly 
spherical  bubbles  are  common  in  nature  and  convenient  to  work  with  in 
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experiments.  As  a  consequence  of  their  symmetry,  such  bubbles 

display  a  scattering  enhancement,  known  as  the  glory,  in  the  forward 

and  backward  directions;  critical- angle  effects  like  those  to  be 

described  may  be  observed  from  bubbles  of  other  shapes  as  well. 

Three  models  of  the  scattering  are  of  interest  for  comparison  with 

observations:  the  Mie  theory,  the  geometrical-optics  approach,  and 

physical-optics  models.  These  will  be  discussed  in  Section  1.2. 

The  optical  properties  of  bubbles  are  of  Interest  in  a 

variety  of  fields.  In  oceanography,  for  example,  the  vast  numbers  of 

gas  bubbles  which  inhabit  the  upper  layers  of  the  ocean  are  seen  as 

3  4 

important  components  of  the  ecosystem,  '  as  well  as  being  an 

influence  on  the  color  of  the  sea^  and  on  its  sound  transmission 
6 

properties.  Various  probes  have  been  used  to  measure  the  oceanic 

6-8 

bubble  size  distribution,  but  considerable  disparity  exists 
between  optical  and  acoustical  findings;  it  appears  that  more 

9 

sensitive  optical  techniques  are  needed.  In  the  science  of  nuclear 
particle  detection,  the  analysis  of  tracks  in  bubble  chambers  may  be 
aided  by  more  sophisticated  models  of  the  scattering  from 
bubbles.^  ^  In  hydrodynamics  research,  bubbles  have  been  used  as 

13 

scatterers  for  laser-Doppler  velocimetry,  and  as  tracers  in 

14 

photographic  flow-visualization  studies.  The  work  presented  in 
this  thesis  may  be  of  value  in  any  application  where  bubbles  are  to 
be  detected,  sized,  or  discriminated  from  other  types  of  particles. 
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1«2  Scattering  Modal* 

A)  The  Mie  Scattering  Theory 

A  mil-known  exact  solution  exists  for  the  scattering  pattern 
'of  a  homogeneous  dielectric  sphere  under  plane-wave  illumination. 
Published  by  Mie^  in  1908,  the  solution  is  presented  in  a  modern 
form  by  many  authors  (e.g.  Refs.  16-18).  The  Mie  theory  has  found 
widespread  application  to  scattering  problems  since  it  places  no 
formal  restrictions  on  the  size  or  the  refractive  index  of  the 
sphere.  Some  details  of  the  Mie  solution  which  are  relevant  to  the 
problem  at  hand  will  be  reviewed  here. 

Let  the  sphere  have  a  radius  a  and  relative  refractive 
index  m,  and  let  the  incident  light  be  of  wavelength  Xq  in  the  outer 
medium.  The  dimensionless  size  parameter  commonly  used  is  ka  - 
2sa/Ao,  where  k  is  the  wavenumber  of  the  light;  ka  »  1  for  all 
the  bubbles  observed  in  this  research.  The  Mie  theory  gives  the 
amplitude  and  phase  of  the  scattered  electromagnetic  field  inside  and 
outside  the  sphere,  but  in  the  present  work  only  the  field  at 
distances  R  »  a  from  the  center  of  the  sphere  is  of  interest. 

The  scattered  light  from  a  sphere  can  be  completely 
characterized  by  the  complex  amplitude  functions  S^(<J>)  and  Sj C  d>  ) , 
where  $  is  the  scattering  angle  measured  with  respect  to  the  incident 
wave's  direction.  The  subscripts  1  and  2,  respectively,  denote  the 
polarization  components  perpendicular  and  parallel  to  the  scattering 
plane;  the  plane  of  the  scattering  is  defined  by  the  directions  of 
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the  incident  and  scattered  light.  At  a  distance  R  the  intensity  of 
the  scattered  radiation  having  polarization  j  ■  1  or  2  is 

IjCR,*)  -  I0j|Sj|2(kRr2  U-D 

where  Iqj  is  the  intensity  of  the  j-polarized  component  of  the 
incident  light.  Most  often  in  this  dissertation  the  scattered 
intensity  will  be  normalized  to  the  intensity  I_  from  a  perfectly 
reflecting  sphere  of  radius  a, 

IR  -  Ioj(a/2R)2  2  1,  (1-2) 

in  which  case  the  scattered  intensity  Ij  is  given  by 

Ij(K,*)  -  | S j | 2  (2/ka)2.  (1.3) 


The  Mie  theory  gives  the  scattering  amplitudes  Sj($)  in  the  forms 


s1(<t») 


2n  +  1 


n  *  1 


(a"  +  b  T  ), 


n(n  +  1)  'an  T  Vn 


¥*>-  I  ItsVtt  <Vn  ♦  *.v- 

n  *  i 

The  angle  dependence  in  Eq.  (1.4)  is  contained  in  the  functions 


(1.4) 


*n  ■  P*(cos  #)/sin  $ 
Tn  -  d[P*(coa  $)]/  d$. 


(1.5) 


using  the  associated  Legendre  polynomial 


The 


P^Ccos  ♦). 

coefficients  a  and  b  in  Eq.  (1.4)  are  functions  of  the 
n  n 

physical  parameters  of  the  problem.  Letting  x  ■  ka  and  y  -  mka, 
these  coefficients  are  expressed  by 

ra(y)  fn(*)  -  ■Vy)fi(*) 
a“  “  ’n(y)5n(5t)  "  BVy)5i(x) 

(1.6) 

nr(y)4'n(x)  -  fn(y)^(x) 

bfl  "  n»;(y)Cn(x)  -  Vy)5n(x) 

using  the  Ricatti-Bessel  functions  and  their  derivatives 


¥„(*)  •  (iw*)^  Jn+^(z), 
r(z)  -  df  (z)/dz, 

n  n 

Cn(z)  -  airz)^Hn+^(z), 

?„(z)  -  d?n(z)/dz, 


where  J^^z)  is  the  Bessel  function  of  the  first  kind  and  Hq+|  (2) 
is  the  Hankel  function  of  the  first  kind. 

The  solution  outlined  above  is  correct  for  spheres  of 
arbitrary  size.  There  is,  however,  a  practical  limitation  on 


Mie- theory  computations  that  arises  in  connection  with  the  size 


parameter  ka.  The  expressions  (1.4)  for  the  Sj($)  converge  after 
n  a  ka  terms  of  the  series  so  that  for  a  single  angle  $  the 
computing  time  Increases  as  ka.  Moreover,  when  scattering  results 
are  desired  over  a  fixed  range  of  angles,  the  angular  resolution  of 
the  computations  may  also  have  to  be  increased  roughly  in  proportion 
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19 

to  lea  to  allow  floe  details  to  be  resolved .  The  steps 

(1.5)-(1.7)  required  to  generate  each  term  of  Eq.  (1.4)  are  somewhat 
20 

time-consuming ,  hence  for  large  values  of  the  size  parameter  Mle 

calculations  Involve  a  considerable  Investment  In  computing  time. 

Though  they  cannot  dispel  this  basic  difficulty,  efficient  algorithms 

19 

such  as  those  of  Wlscombe  are  welcomed  when  the  Hie  theory  is  used 
to  study  angular  scattering  features  or  the  effects  of  changes  in  the 
parameters  m  or  ka. 

It  should  be  noted  that  the  assumption  of  an  incident  plane 
wave  is  not  regarded  as  an  inherent  restriction  on  Mle  results.  This 
assumption  is  a  common  one  in  scattering  theories  since  any 
electromagnetic  wave  can  be  decomposed  into  a  linear  combination  of 
plane-wave  components  by  Fourier  analysis.  Similarly,  the  scattered 
field  arising  from  a  non-planar  incident  wave  can  be  regarded  as  a 
superposition  of  the  scattered  fields  from  incident  plane  waves.  But 
where  such  a  solution  is  called  for,  the  calculation  of  Mie  results 
to  be  superposed  means  yet  another  considerable  increase  in  computing 
time. 

Exact  solutions  to  the  scattering  problem  are  also  available 
for  particles  which  are  not  homogeneous  spheres;  a  number  of  these 
are  discussed  by  Bohren  and  Huffman  (Ref.  18,  Chap.  8).  Different 
shapes  or  other  properties  of  the  particle  introduce  additional 
parameters  into  the  theory,  of  course,  and  make  the  computations  much 
more  arduous  than  for  the  Mle  equations.  4a  previously  mentioned, 
the  bubbles  observed  in  the  research  presented  here  were  very  nearly 


spherical;  hence  no  "improved1*  exact  models  of  the  scattering  were 
considered  necessary. 

It  is  possible  to  use  exact  scattering  solutions  to  study  the 
effects  of  changes  in  the  size,  shape,  or  refractive  index  of  the 
scatterer ,  and  of  variations  in  the  form  of  the  incident  wavefront . 
But  approximate  models  which  will  account  for  the  resulting  trends 
become  rather  attractive  in  view  of  the  computational  difficulties 
mentioned  above  for  exact  solutions,  especially  when  the  scattering 
particle  is  large.  At  the  same  time,  however,  the  exact  solution 
allows  the  validity  of  such  approximations  to  be  tested  without 
introducing  the  uncertainties  that  go  with  experiments.  The  Mie 
theory  is  a  valuable  aid  in  developing  models  and  understanding  the 
physics  of  light  scattering  from  bubbles. 

(B)  Geometrical  Optics 

A  useful  approximate  solution  to  a  given  scattering  problem 
can  often  be  developed  based  on  simple  geometrical  optics.  In  this 
approach  the  incident  light  is  treated  as  consisting  of  rays  whose 
scattering  angles  and  intensities  follow  from  the  laws  of  reflection 
and  refraction  applied  at  the  surface  of  the  scattering  object.  The 
total  intensity  in  any  direction  is  then  obtained  by  summing  the 
intensities  of  individual  scattered  rays.  This  type  of  solution 
neglects  interference  and  diffraction  phenomena,  usually  with  the 
assumption  that  these  effects  are  small  when  the  scatterer  is  large 
compared  to  the  wavelength  of  the  light.  This  assumption  is  not 
always  correct,  but  the  simple  model  may  still  give  some  acceptably 


accurate  results  and  some  useful  insights.  As  a  well-known  example, 

the  existence  of  the  rainbow  and  several  of  its  main  features  can  be 
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understood  from  the  geometrical  optics  of  sunlight  on  a  water  drop. 

A  complete  theory  of  the  rainbow  must  include  diffraction,  however, 

since  the  naive  model  incorrectly  ascribes  to  it  an  infinite 

intensity.  For  most  other  angles  the  range  of  wavelengths  in 

sunlight  will  mask  the  effects  of  Interference  so  that  geometrical 

optics  can  approximate  the  average  scattered  intensity.  But  where 

diffraction  is  prominent.  Interference  effects  will  also  appear;  for 

water  drops  this  includes  the  vicinity  of  the  rainbow  and  also  the 

forward  and  backward  directions.^  The  rainbow  exists  only  for 

scatterers  which  are  droplike  (  m  >  1  ) .  For  bubbles  (  m  <  1  )  total 

reflection  of  certain  incident  rays  will  occur,  giving  rise  to  new 

diffraction  effects  in  what  is  called  the  critical-angle  scattering 

region.  The  forward  and  backward  directions  also  exhibit  diffraction 

effects  for  bubbles;  as  is  the  case  for  drops,  the  scattterlng  in 

these  regions  is  extraordinarily  bright  and  is  frequently  referred  to 

as  the  "glory”.  A  geometrical-optics  approximation  was  given  by 
22 

Davis  for  the  scattering  from  a  large  air  bubble  in  water.  Some 
aspects  of  this  kind  of  model  will  be  reviewed  here  with  emphasis  on 
the  regions  of  diffraction,  where  the  need  for  a  more  sophisticated 
approach  becomes  apparent. 

The  goal  of  the  geometrical  optics  model  as  defined  above  is 
to  approximate  the  scattered  lntemaity  as  a  function  of  the 
scattering  angle  ♦  .  First,  it  most  he  determined  which  incident  rays 
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are  scattered  to  a  given  4>.  Next,  the  intensity  associated  with  each 
of  these  rays  oust  be  found.  The  bubbles's  cross-section  in  the 
scattering  plane  is  assumed  to  be  circular  with  radius  a,  and  its 
refractive  index  is  m.  Figure  1.1  illustrates  several  rays  directed 
to  the  same  scattering  angle.  It  is  helpful  to  denote  each  ray  by  a 
pair  of  parameters  (p,  £,  ) ,  where  p  «  the  numher  of  chords  the  ray 
makes  inside  the  bubble  and  1  -  the  number  of  times  the  ray  crosses 
the  optic  axis  (including  crossings  which  may  occur  after  the  ray 
exits  the  bubble).  The  scattering  angle  for  a  (p,£)  ray  is 

<J»  -  2(-l)A[ppp  -  9p  +  i»(i  -  P  +  Cz)],  (1.8) 

where  E^  ■  1  +  (-1)  ]  and  0p  and  0p  are  the  positive  local 
angles  of  incidence  and  refraction,  respectively.  For  any  <J> ,  an 
infinite  number  of  (p,Jl)  rays  will  exist,  subject  to  Snell's  law, 


sin  6  -  m  sin  p  . 

P  P 


(1.9) 


For  a  given  value  of  $,  the  incidence  angle  9  of  a  (p,£)  ray  will 

P 

generally  be  found  numerically  from  Eqs.  (1.8)  and  (1.9).  However, 

exact  trigonometric  solutions  for  0p($)  are  possible  in  some  cases 

(see  Appendix  B  of  Chapter  3). 

The  intensity  of  a  (p,&)  ray  is  reduced  at  each  Internal 
2 

reflection  by  a  factor  r^  , 


where 
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Figure  1.1.  Rays  in  the  scattering  plane  emerging  at  the  scattering 
angle  $  ■  50°.  Rays  are  specified  by  the  parameters  p  -  the  number 
of  chords  within  the  bubble ,  and  l  •  the  number  of  crossings  of  the 
optic  axis. 
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3in(Qp  -  P0) 

1  sin(0p  +  pp) 


tan(0  -  p  ) 

P _ 

^  tan(0p  +  Pp) 


(1.10) 


are  Che  Fresnel  reflection  coefficients  for  the  light  polarized 

perpendicular  (j  *1)  and  parallel  (j  ■  2)  to  the  plane  of  the 

scattering .  Transmissions  Into  and  out  of  the  bubble  Introduce  a 
2  2 

factor  (1  -  r  )  into  the  scattered  intensity.  Also,  a  geometric 
divergence  factor  must  be  included  to  account  for  the  spreading 

of  the  incident  energy  over  a  different  ares'  after  the  scattering. 
The  geometry  of  the  problem  leads  to  the  expression 


D  -  ^  sin  20  |d0  /dd»|  /  sin  $  (1.11) 

r  r  r 

(see  Ref.  22,  Sec.IV-A).  From  Eqs.  (1.8)  and  (1.9)  above. 


|d<D/d0p|  -  2|pr  -  1| 


(1.12) 


where  T*dp/d0  *  tan  p  /  tan  0  ,  hence 
P  P  P  P 


D  ■  £  sin  20  I PT  -  l|”Vsinp. 

p  p<«- 


(1.13) 


By  this  analysis,  the  intensity  of  a  (p,£)  ray  is  found  to  be 


(1.14) 
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normalized  to  the  intensity  from  a  perfectly  reflecting  sphere  as  in 
Eq.  (1.3).  The  total  scattered  intensity  at  $  according  to  simple 
geometrical  optics  is  then 


I,(*)  -  I 

J  <p,4)  PJ 


(1.15) 


The  summation  in  Eq.  (1.15)  includes  all  possible  (p,2  )  rays 

scattered  to  <J> ,  but  ordinarily  only  those  rays  having  small  values  of 

p  are  necessary,  because  of  the  factor  rj^P  ^  in  Eq.  (1.14). 

The  contribution  from  a  ray  with  many  internal  reflections  becomes 

insignificant,  except  when  becomes  very  large,  as  noted  below. 

There  are  certain  angles  in  the  scattering  from  bubbles  which 

require  special  attention,  as  further  examination  of  will  show. 

First,  for  axial  rays  (incident  along  the  optic  axis),  0  *  0  and  *  » 

P 

0  or  180°,  so  that  Eq.  (1.13)  is  indeterminate.  But  its  limiting 

-1  -2 

form  in  this  case  becomes  Dp  »[2(pm  -1)]  ;  axial  rays  present 

no  particular  difficulties  if  this  form  for  is  used.  There 

will,  however,  also  be  rays  scattered  to  angles  <t>  ■  0  or  180°  but 

with  9  *  0.  In  such  cases  D  becomes  infinite.  These  are  the 

P  P 

glory  rays,  scattered  parallel  to,  but  not  coincident  with,  the  optic 
axis.  Their  relative  intensities  may  be  quite  large  but  they  are  not 
infinite;  glory  scattering  cannot  be  properly  dealt  with  by  simple 


geometrical  optics. 


Another  exceptional  angle  is  evident  in  graphs  of  the 
Intensity  distribution  obtained  in  Eq.  (1.15)  (e.g.  see  Fig.  18  of 
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Ref.  22  or  Fig.  5  of  Ref.  12).  A  cusp  occurs  in  the  intensity 

pattern  at  the  critical  scattering  angle  <t>c  ■  130°  -  2  8^,  where  8^  * 

arcsin  m  is  the  critical  angle  of  incidence.  From  Eq.  (1.9)  the 

refraction  angle  becomes  90°  when  8  *  8  ,  and  so  T  »°°.  Then  the 

only  ray  having  a  nonzero  value  for  D^  is  the  externally  reflected 

(0,0)  ray. Rays  incident  at  angles  >_  9^  are  totally  reflected  at  the 

bubble's  surface,  contributing  to  an  enhanced  intensity  in  the  region 

from  (J>  ■  0  to  <t>  .  The  transition  to  total  reflection  also  introduces 
c 

a  phase  shift  into  the  scattered  ray,  since  the  reflection 

coefficients  in  Eq.  (1.10)  become  imaginary  when  8^  exceeds  8^. 

Simple  geometrical  optics  is  inadequate  for  describing  this 

transition  region.  The  angular  derivative  of  the  intensity,  dl/d<t> 

is  predicted  to  become  infinite  at  <&  ,  which  is  physically  untenable, 

and  Eq.  (1.15)  does  not  succeed  in  predicting  even  the  average 

intensity  for  angles  in  the  near  vicinity  of  <t>c. 

The  geometrical-optics  approach  gives  helpful  insights  about 

the  scattering  properties  of  bubbles,  and  also  serves  to  identify  the 

angles  where  a  more  sophisticated  model  is  necessary.  The  existence 

of  the  critical- angle  scattering  region  is  a  consequence  of  the 

refractive  index  m  being  less  than  unity  for  bubbles.  Details  of 

the  phenomena  which  appear  in  this  region  were,  until  recently,  not 

well-explored  and  will  be  the  topic  of  Chapter  2.  Glory  scattering 

is  found  to  be  a  result  of  the  spherical  symmetry  of  the  bubble. 

16  23 

Analogous  effects  exist  in  the  scattering  from  spherical  drops,  ’ 

’4 

and  also  in  particle  scattering  by  central  forces."  Chapter  3  will 
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present  models  and  measurements  of  the  phenomena  which  accompany  the 
forward  and  backward  glory  of  bubbles. 

(C)  Physical  Optics 

The  physical-optics  approximations  used  in  this  research  are 
described  in  detail  in  the  following  chapters  and  their  references. 
Here  it  Is  worthwhile  to  note  some  general  characteristics  of  the 
method  for  the  purpose  of  comparing  it  with  the  two  models  described 
above . 

Unlike  simple  geometrical  optics,  the  physical-optics 
approach  recognizes  the  wave  nature  of  light  and  deals  with 
diffraction  and  interference  effects  in  the  scattering.  But  unlike 
the  Mle  theory,  full  solutions  to  Maxwell’s  equations  are  not 
determined.  Instead,  the  physical-optics  model  approximates  certain 
scattered  wavefronts  near  the  bubble  and  uses  diffraction  theory  to 
find  the  resulting  pattern  of  light  far  away.  Much  of  the  physical 
insight  given  by  the  geometrical  approach  is  retained,  however,  in 
formulating  characteristics  of  the  scattered  waves.  The  curvature 
and  relative  phases  of  wavefronts  are  determined  geometrically,  while 
their  amplitudes  are  derived  using  the  Fresnel  coefficients  and 
appropriate  divergence  factors.  Sy  considering  only  the  most 
significant  scattered  waves,  the  physical-optics  approach  can  give  a 
fair  approximation  to  features  in  the  scattering  while  avoiding  the 
computational  complexities  of  the  Mia  theory.  The  physical-optics 
approximations  presented  here  are  success  '■  In  removing  the 
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unphysical  divergences  mentioned  above,  and  they  provide  Insights 
into  the  scattering  process. 
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CHAPTER  2 


CRITICAL-ANGLE  SCATTERING 


2.1.  Introduction 

Observations  of  light  scattered  near  the  critical  angle  by 

air  bubbles  in  water  reveal  an  interesting  oscillatory  structure  in 

the  far-zone  intensity.  Monochromatic  illumination  of  a  single 

bubble  gives  rise  to  coarse  undulations  of  intensity,  reminiscent  of 

an  edge  diffraction  pattern,  along  with  a  superposed  fine  structure.* 

With  white  light  the  fine  structure  is  lost,  while  the  coarse 

structure  is  manifested  as  a  series  of  colored  bands;  this  has  been 

2  3 

observed  for  a  cylindrical  bubble  in  glass.  Pulfrich  (in  1888) 
reported  rainbow-like  colors  from  a  sunlit  cloud  of  small  bubbles 
rising  in  water;  his  observations  were  made  with  the  unaided  eye  in 
Che  critical  scattering  region.  The  present  chapter  gives  Che  first 
detailed  experimental  investigation  of  structures  in  the  scattering 
from  rising  bubbles.  Features  observed  in  the  scattering  of  laser 
light  by  single  bubbles  will  be  compared  with  Che  predictions  of 
models. 

Figure  2.1(a)  illustrates  three  models  of  scattered  intensity 

computed  near  the  critical  angle.  The  solid  line  in  this  figure 

4 

represents  the  Mie  theory,  which  gives  the  exact  scattering  pattern 
for  a  dielectric  sphere  under  plane-wave  illumination;  the  coarse  and 
fine  structures  are  plainly  visible.  The  dashed  line  is  from  a 
physical-optics  model** ^  which  has  proven  successful  in  approximating 
the  coarse  structure  in  Mie  results  over  a  range  of  bubble  sizes.*’ ^ 


2. 


Figure  2.1.  Normalized  scattered  intensity  from  a  bubble 
with  ka  ■  1633  and  the  electric  field  parallel  to  the 
scattering  plane,  (a)  Three  models:  the  solid  curve  is 
the  Mie  theory,  the  dashed  curve  is  a  physical-optics 
approximation,  and  the  dotted  curve  is  from  simple 
geometric  optics,  (b)  The  solid  curve  is  experimental 
data  taken  from  the  photograph  in  Fig.  2.6(a),  and  the 
dashed  curve  is  the  physical-optics  model. 
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Ic  applies  far-zone  diffraction  theory  to  certain  scattered  rays 
whose  amplitudes  adjacent  to  the  bubble  are  computed  using 
plane-surface  reflection  coefficients.  The  dotted  line  comes  from 

g 

the  naive  geometric-optics  model  of  Davis;  It  neglects  the  effects 
of  diffraction  and  interference,  resulting  in  a  divergent  derivative 
of  intensity*  at  the  critical  angle  (  2  82.65°  for  this  example)  and 
the  loss  of  all  oscillatory  structure.  Figure  2.1(b)  shows 
experimental  data  plotted  along  with  the  physical-optics 
approximation.  Details  of  the  experiment  and  normalization  of 
results  will  be  discussed  In  Secs.  2.3  and  2.6,  but  notice  here  how 
the  data  compare  with  the  different  models.  There  is  a  strong 
similarity  to  the  Hie  result  in  the  appearance  of  the  coarse  and  fine 
structures.  The  average  intensity  near  the  critical  angle  does  not 
exhibit  the  divergent  derivative  predicted  by  simple  geometric 
optics,  but  Instead  decays  gradually  as  in  the  physical-optics  model 
and  Mle  theory. 

In  this  work  attention  has  been  limited  to  bubbles  that  are 
much  larger  than  the  wavelength  of  the  Incident  light.  While  the  Mie 
theory  is  certain  to  give  the  most  complete  description  of  the 
scattering  for  a  spherical  bubble,  it  gives  no  direct  insight  into 
the  origin  of  observed  features  in  the  scattering,  nor  into  the 
effects  of  nonsphericity,  a  condition  to  be  expected  in  dealing  with 
real  bubbles  of  large  size.  The  physical-optics  approximation 
elucidates  the  roles  of  certain  rays  in  producing  the  coarse 
structure,  and  a  similar  approach  is  helpful  in  understanding  the 
fine  structure  as  well.  It  may  also  be  extended  into  the  realm  of 
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nonspherlcal  bodies  (e  physical-optics  approximation  has  been  applied 

a 

to  large,  spheroidal  drop-like  scatterers  ). 

Consider  a  plane  wave  incident  upon  a  sphere  of  refractive 
index  m  -  n^/nQ,  where  the  indices  of  the  inner  and  outer  media, 
n^  and  nQ,  are  real,  and  <  nQ.  The  size  parameter  commonly 

used  is  ka  •  2ira/XQ,  where  a  is  the  sphere  radius  and  k  is  the 
wavenumber  for  the  incident  light.  Note  that  XQ,  the  wavelength  in 
the  surrounding  medium  (water),  is  related  to  the  vacuum  wavelength 
by  XQ  ■  Let  denote  the  scattering  angle,  measured 

between  the  direction  of  the  incident  wave  propagation  and  that  of  a 
scattered  ray  leaving  the  bubble.  Figure  2.2  illustrates,  for  m  » 
3/4,  several  different  rays  scattered  to  the  same  angle  $  ■  50 0  . 
These  rays  may  be  characterized  by  the  parameters  p  -  the  number  of 
chords  the  ray  makes  within  the  bubble,  and  1  -  the  number  of  times 
the  ray  crosses  the  optic  axis.  In  this  chapter  rays  will  often  be 
specified  by  their  parameters  (p,  l  ) .  The  scattering  angle  for  a 
(p,Z)  ray  will  be  given  by 

*  -  2(-l)*[ppp  -  8p  +  *ir(i  -  p  +  e4)],  (2.1) 

o 

where  e«  •  (1  +  (-1)  )/2,  and  0  and  p  are  the  positive  local  angles 

*  p  p 

of  incidence  and  refraction,  respectively.  All  rays  satisfy  Snell's 
law. 


sin  •  m  sin  p 
P  P 


(2.2) 
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Figure  2.2.  Rays  in  the  scattering  plane  emerging  at  the  scattering 
angle  $  ■  50°.  Rays  are  specified  by  the  parameters  p  -  the  number 
of  chords  within  the  bubble,  and  l  -  the  number  of  crossings  of  the 
optic  axis. 
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At  tha  critical  angle  of  incidence  0c  ■  arcsin  a  the  refraction 

angle  becoaes  90*.  For  the  equivalent  plane-interface  problem,  the 

reflection  is  total  when  the  Incidence  angle  exceeds  3  .  In  the 

c 

physical-optics  approximation  of  the  scattering ,  total  reflection 
occurs  for  0^  >  9£  provided  ka  »  1  so  that  tunneling  through  the 
bubble  is  negligible. 

Because  of  this  onset  of  total  reflection,  there  is  a 

critical  scattering  angle  *  x  -  29^  above  which  no  refracted  1*0 

ray  may  exist.  But  along  with  the  externally  reflected  (0,0)  ray 

there  will  be  many  rays  having  1  >  0  directed  to  angles  >  4>c .  Of 

these,  the  one  contributing  the  greatest  intensity  according  to  ray 

0 

optics  will  be  the  (2,1)  ray,  which  suffers  only  one  internal 
reflection.  The  physical-optics  model  of  the  coarse  structure  makes 
use  of  the  (0,0)  and  (1,0)  rays  to  compute  the  far-zone  intensity  in 
the  vicinity  of  <l»c  (see  Appendix  A  for  a  summary  of  the  model).  In 
the  present  chapter  the  coarse  structure  observed  from  real  bubbles 
will  be  compared  with  the  predictions  of  this  model.  It  will  also  be 
shown  that  the  interference  of  the  (0,0)  and  (2,1)  rays  is  able  to 
approximate  the  observed  fine- structure  period  in  the  critical 
region.  In  addition,  the  effect  of  the  (3,1)  ray  on  the  fine 
structure  will  be  considered. 

Critical-angle  scattering  phenomena  will  be  present  for  any 
scatterer  with  relative  refractive  index  m  <  1.  An  understanding  of 
these  phenomena  should  be  advantageous  in  the  design  of  optical 
instruments  to  detect  bubbles,  size  them,  or  discriminate  them  from 
particulate .  Applications  may  include  measurements  of  microbubble 
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populations  In  the  sea,^  the  photography  of  hydrogen  tracer  bubbles 

for  hydrodynamics  research,  laser-Doppler  velocimetry  of  bubbles 

11  12 
carried  by  liquids,  and  the  detection  of  bubbles  in  glass.  Many 

of  these  results  will  also  pertain  to  the  analogous  problem  in 

13 

acoustical  scattering . 

2.2.  Calculation  of  the  Fine  Structure  Period 

When  a  plane  wave  is  incident  upon  a  bubble,  the  scattered 
rays  will  have  curved  wavefronts,  as  if  they  were  emerging  from 
source  points  in  the  scattering  plane.  To  locate  the  virtual  source 
point,  consider  a  (p,£)  ray  scattered  to  an  angle  <!> ,  and  a  similar 
ray  in  the  same  plane  scattered  to  d>* ;  Che  backward  extrapolations  of 
these  two  rays  will  intersect.  The  virtual  source  point  Fp  is 
defined  by  this  intersection  in  the  limit  as  p'  P  .  Figure  2.3 
shows  the  virtual  sources  Fq  and  for  the  (0,0)  and  (2,1)  rays. 

Let  <Xp  be  the  distance  from  Fp  to  an  exit  plane  which  is  tangent  to 
the  bubble  and  normal  to  the  direction  of  scattering .  It  is  shown  in 
Appendix  B  chat 


ap  -  a[  1  +  j(pT  -  I)_1cos0p],  (2.3) 

with  T  ■  tan  pp/tan  9p.  Because  9  and  p  change  with  p,  the  location 
of  each  Fp  varies  with  the  scattering  angle.  To  find  the  angular 
period  of  interference  Ap  for  two  scattered  rays  while  Including  the 
movement  of  their  virtual  sources,  it  is  helpful  to  consider  their 
relative  phases. 
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Figure  2.3.  Rays  used  to  model  the  fine-structure  interference 
period.  Points  F0  and  F2  are  the  virtual  source  locations, 
respectively,  of  the  rays  labeled  0  and  2  in  the  limit  as  the 
dashed  O’  and  2’  rays  approach  them,  and  Op  is  the  distance  from 
F  to  the  exit  plane. 


Letting  the  phase  of  the  incident  plane  wave  be  zero  at  the 
entrance  plane  tangent  to  the  bubble  and  normal  to  the  wave's 
direction,  the  phase  at  the  exit  plane  may  be  specified  for  any  (p,£) 
ray.  First,  there  is  a  phase  delay  due  to  propagation  through  the 
media  from  the  entrance  plane  to  the  exit  plane,  given  by 

np(ep'pp)  *  2  ka  (1  -  cos  9p  +  pm  cos  pp),  (2.4) 

where  0^  and  are  the  incidence  and  refraction  angles  of  the  (p,&  ) 

ray.  There  will  also  be  a  phase  advance  5^  associated  with  total 

reflection  of  the  (0,0)  ray;  for  >  9C»  the  Fresnel  reflection 

14 

coefficients  give 

tan  »  m^*  -^(sin^Q  -  m^)^/cos  9q,  (2.5) 

where  the  subscript  j  is  assigned  the  value  l  when  the  incident 
electric  vector  is  entirely  perpendicular  to  the  scattering  plane  and 
the  value  2  when  entirely  parallel.  There  also  occurs  a  phase 
advance  of  ir/2  whenever  a  ray  crosses  a  focal  point. The  two  types 
of  focal  points  present  in  the  scattering  from  bubbles  can  be  seen  in 
Fig.  2.3.  One  focus  occurs  where  adjacent  rays  (such  as  the  nearby 
two-chord  rays  labeled  2  and  2')  cross  in  the  scattering  plane;  a  ray 
with  p  >  0  chords  experiences  p  -  1  such  crossings  with  adjacent 
rays.  The  point  where  a  ray  crosses  the  optic  axis  is  also  a  focus 
since  similar  rays  in  other  scattering  planes  cross  the  axis  at  the 
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same  point.  Thus,  a  phase  advance  of  (p  -  1  +2,)  n/2  is  incurred  by 
a  (P»2)  ray  (with  p  >  0)  due  to  passage  through  focal  points. 

For  the  (0,0)  and  (2,1).  rays  scattered  to  angle  <J>,  the  phase 
difference  at  the  exit  plane  will  be  [^(^.f^)  -  it]  -  [HqOq.Pq)  - 
6  (0O)].  For  a  similar  pair  of  rays  scattered  to  <J>'  -  $  -  A<j>,  the 
phase  difference  will  be  identical  in  form  but  with  all  angles 
primed.  Now,  if  one  period  of  the  far-zone  interference  pattern  of 
the  (0,0)  and  (2,1)  ray3  occurs  from  <J>  to  ,  the  phase  difference 
must  change  by  2ir.  The  condition  for  one  interference  period  in  the 
angular  spread  A<J>  becomes 

2t r  -  2  ka{  [cos  0,(1  -  cosA9,)-  sin  9,  sin  A0,  ](2J  -  1) 

Z  (2.6) 

+  cos  0Q(  1  -  cos  A9q)  -  sin  9q  sin  A0Q}  +  5j(0Q)  -  6 ^ (0^) , 

where  J  «  tan(  tan(  A02/2),  A0p  -  8p  -  0  '  ,  and  A^  -  02  - 

02'«  To  express  this  condition  in  terms  of  A4> »  the  expression  for  5^ 
may  be  approximated  in  the  vicinity  of  0^  by 

6j(0o)  »  ■2(1_j)(8tan8c)i(90  -  0(;)i  (2.7) 

for  0Q  >_  0c,  and  by  5  -  0  for  0Q  <  0^.  From  Eq.  (2.1),  A<t>  -  2  A&,  - 
4Ap2  -  -2A0Q.  Making  these  substitutions  in  Eq.  (2.6)  and  retaining 
terms  down  to  order  A#  give  a  quadratic  equation  in  ( A<fc  )  the 
resulting  fine-structure  angular  period  (in  radians)  is 

A4»  -  Xq/B2  -  0[(Xo/B2)3/2], 


(2.8) 


where  ■  a(sin  0^  +  sin  is  the  sum  of  the  impact  parameters 

of  the  two  rays  and  also  their  lateral  spacing  in  the  exit  plane. 

To  a  first  approximation  the  fine-structure  period  is 

equivalent  to  the  interference  period  of  two  point  sources  separated 

3 /2 

by  a  distance  B^.  The  0{(  J^/B^)  J  term  comes  from  including  the 

reflection  phase  shift  of  the  (0,0)  ray  when  0g  >  0C>  while  terms 

attributable  to  movement  of  the  virtual  sources  of  the  scattered  rays 

2 

are  of  the  order  and  3tna^er  •  ln  Sec.  2.5  of  this  chapter 

the  first-term  approximation  of  Eq.  (2.8), 

a  X  /B0,  (2.9) 

0  z 

is  shown  to  give  a  fair  description  of  Che  fine-structure  period 

observed  from  real  bubbles .  This  approximation  was  previously 

7  13 

obtained  using  elementary  methods.  ’ 

2 .3 .  Experimental  Methods 

The  experiments  that  are  described  in  this  section  provided 
absolute-angle,  relative  intensity  measurements  of  the  far-zone 
scattering  from  single  bubbles  in  the  critical  angle  region.  To 
minimize  the  effects  of  bubble  nonsphericity,  the  scattering  plane 
was  chosen  to  be  horizontal.  The  chamber  in  which  the  scattering 
took  place  was  fashioned  from  an  aluminum  block  with  portholes,  and 
was  filled  with  distilled  water.  The  incident  light,  from  a  He-Ne 
laser,  entered  the  chamber  through  a  front  window,  and  the  scattering 
was  photographed  by  a  camera  through  a  side  window.  The  camera  lens 
was  focused  at  infinity  so  that  the  far-zone  scattering  pattern  was 
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recorded.  Bubbles  were  injected  by  a  syringe  and  a  steel  needle 
which  entered  through  the  floor  of  the  chamber.  The  scattering  was 
photographed  either  with  the  bubble  in  a  pendant  form  (attached  to 
the  tip  of  the  needle),  or  else  rising  freely  through  the  water. 
Bubble  diameters  were  measured  using  a  microscope  situated  in  air 
above  the  top  opening  of  the  scattering  chamber.  Rising  bubbles  were 
trapped  by  a  glass  slide  across  this  opening  so  their  sizes  could  be 
measured.  From  microdensitometer  scans  of  negatives,  relative 
intensity  profiles  could  be  derived  for  the  scattering  by  bubbles. 

Since  the  camera  was  in  air,  refraction  at  the  viewing  window 
had  to  be  taken  into  account  in  determining  the  absolute  scattering 
angles  represented  in  the  photographs.  Figure  2.4  shows  a  diagram  of 
the  apparatus  as  it  was  used  to  perform  the  angle  calibration.  A 
flat  glass  reflector  was  suspended  vertically  through  the  top  opening 
of  the  water-filled  chamber.  By  turning  the  reflector  about  a 
vertical  axis,  the  incident  laser  beam  could  be  directed  to  angles  in 
the  horizontal  plane.  A  goniometer  attached  to  the  axis  of  this 
rotatable  reflector  allowed  its  orientation  to  be  measured.  To 
relate  the  goniometer  readings  to  scattering  angles,  the  setting 
which  gave  direct  backscattering  (  <t>  -  180s)  was  established  using  the 
optical  system  shown  in  Fig.  2.4.  The  backscattering  setting  was 
determined  by  focusing  light  returned  by  the  rotatable  reflector  to 
the  same  spot  on  the  screen  as  light  from  the  beamsplitter  and 
retroreflecting  (corner-cube)  prism.  The  beam  could  then  be  directed 
into  the  camera  at  known  scattering  angles.  A  multiple-exposure 
photograph  of  the  beam  at  several  angles  provided  data  relating 


Figure  2.4(a)  Top-view  diagram  of  apparatus  set  up  for  angle 
calibration.  To  observe  bubble  scattering  the  rotatable 
reflector  was  removed  and  bubbles  were  injected  by  a  needle 
entering  the  bottom  of  the  aluminum  block.  A  microscope  above 
the  top  opening  was  used  to  measure  bubbles.  The  screen  (used 
during  angle  calibration  of  the  rotatable  reflector)  is  at  the 
optical-transform  plane  of  the  lens.  The  systematic  uncertainty 
in  the  angle  calibration  was  <  0.02°. 

Figure  2.4(b)  Photograph  of  the  experimental  apparatus.  The 
numbered  devices  are:  (1)  front  end  of  laser,  (2)  beam  expander, 
(3)  polarization  rotator,  (4)  mirror,  (5)  beam  splitter,  (6) 
retroreflecting  prism,  (7)  lens,  (8)  scattering  chamber,  (9) 
camera,  (10)  mechanism  of  angle-calibration  goniometer  and 
rotatable  reflector. 


33 


points  across  the  camera’s  field  of  view  to  scattering  angles  <^. 
Typically,  about  15  of  these  calibration  points  were  taken  across  a 
negative,  and  each  roll  of  film  contained  one  to  three  such  frames  of 
angle  information. 

The  glass  rotatable  reflector  was  removed  from  the  chamber 

when  scattering  from  bubbles  was  to  be  observed .  The  inner  walls  of 

the  scattering  chamber  had  been  blackened  to  absorb  stray  light.  To 

further  reduce  the  background  intensity,  particulate  larger  chan 

0.2- Um  diameter  were  filtered  from  the  water.  The  incident  beam, 

from  a  15  mW  He-Ne  laser,  had  a  Gaussian  Intensity  profile 

proportional  to  exp(-r^/</*)  where  o*  2.8  mm.  The  beam  was  made  to 

be  highly  collimated  in  the  region  of  the  scattering  using  a 

16 

parallel-plate  shearing  Interferometer .  A  polarization  rotator 
(back-to-back  Fresnel  rhombs)  allowed  the  plane  of  the  beam's 
polarization  to  be  selected;  the  incident  electric  field  was  oriented 
either  (a)  perpendicular  to  the  plane  of  scattering,  corresponding  to 
J  -  1  in  the  preceding  analysis,  or  (b)  parallel  to  the  scattering 
plane,  corresponding  to  j  ■  2.  The  central  region  of  the  incident 
beam  is  expected  to  closely  approximate  a  plane  wave  for  a  scatterer 

whose  radius  is  small  compared  to  the  beam  radius. *  To  place 

pendant  bubbles  in  the  center  of  the  beam,  a  blackened  needle  was 
made  to  stand  vertically  with  its  tip  near  the  beam  axis.  Rising 
bubbles  were  injected  at  the  bottom  of  a  glass  tube  (0.8-cm  i.d.)  and 
rose  into  the  scattering  chamber  where  they  passed  through  the 
central  part  of  the  beam.  A  path  length  of  **  9  cm  helped  dampen 
oscillations  of  the  moving  bubbles.  The  viewing  window,  optically 
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polished  fused  quartz,  was  oriented  so  that  normal  incidence  was  for 
light  scattered  horizontally  at  $3  77*.  For  most  of  the  photography 
the  lens  used  had  an  effective  focal  length  of  85  mm  and  the  viewed 
angular  range  was  limited  to  "*17#;  other  lenses  used  had  effective 
focal  lengths  of  50  and  135  can.  Only  one  lens  was  used  for  each  roll 
of  film,  since  interchanging  lenses  necessitated  a  new  angle 
calibration  as  well. 

Obtaining  properly  exposed  negatives  turned  out  to  be  simpler 
for  rising  bubbles  than  for  pendant  ones.  Each  pendant  bubble  was 
photographed  at  several  shutter  speeds  for  both  polarizations  of  the 
incident  beam.  Best  detail  was  usually  obtained  when  the  camera's 
built-in  light  meter  indicated  a  slight  underexposure;  the  times 
varied  from  ~  1  second  for  small  bubbles  (  -  0.1-mn  radius)  to  1/80 
second  for  large  bubbles  (  -  1-mrn  radius).  With  rising  bubbles  this 
procedure  was  impossible  since  only  a  single  photograph  could  be 
taken  for  each  bubble.  It  was  found,  however,  that  by  simply  holding 
the  camera  shutter  open  as  the  bubble  traversed  the  beam,  acceptable 
exposures  were  obtained  despite  the  wide  range  of  bubble  sizes 
observed.  To  understand  this  fortuitous  result,  consider  a  bubble 
whose  radius,  when  spherical,  is  a.  Its  terminal  velocity  is  found 
to  be  approximately  proportional  ’  to  a  for  a  $  0.7  a,  so 
such  a  bubble  traverses  the  beam  in  a  time  Interval  t  «  a  .  But 

the  scattered  Intensity  If  from  the  bubble,  based  on  geometrical 

8  2 
considerations,  is  also  roughly  proportional  to  a  .  Since  the 

exposure  received  by  the  film  depends  on  the  product  I#t,  it  will  be 

approximately  independent  of  the  bubble  radius  for  small  enough  a 
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values.  For  a  -  0.7  mm  the  terminal  velocity  begins  to  decrease 
2 

from  the  a  dependence,  resulting  in  increased  exposure,  so  an 
appropriate  combination  of  bean  intensity  and  film  sensitivity  may  be 
expected  to  yield  adequate  exposures  up  to  some  bubble  size  limit. 
The  film  used  in  the  experiments  was  Tri-X  (400  ASA).  Rising  bubbles 
were  produced  ranging  in  size  from  ~  0.026-  to  0.81-mm  radius,  and 
acceptable  negatives  were  obtained  throughout  this  range.  For  very 
small  bubbles,  the  background  scattering  often  became  significant 
during  the  long  exposure  times  required .  In  such  cases  photographs 
were  also  taken  with  no  bubbles  present,  to  allow  the  background 

intensity  to  be  subtracted  from  the  scattering  measurements. 

The  absolute  intensity  of  scattered  light  was  not  determined 

in  the  experiments.  Relative  intensity  values  were  important, 

however,  for  comparing  coarse-structure  results  with  models. 

Negatives  were  scanned  with  a  microdensitometer  to  obtain  their 

transmittance  profiles,  and  from  the  transmittance  T  a  value  for 

the  relative  intensity  Ir  was  inferred.  To  make  this  conversion 

from  T  to  possible,  a  response  curve  for  each  film  was 

determined.  A  relationship  exists  between  the  photographic  density 

D  “  -log  T  of  a  negative  and  the  exposure  E  which  it  received, 

where  E  *  I  t,  with  I  the  absolute  intensity  and  t  the  time  of 
a  a 
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exposure.  Film  response  is  usually  shown  by  the  H-D  curve,  where 

0  is  plotted  vs  log  E.  In  the  experiments  several  frames  on  each 

film  were  exposed  to  a  constant  radiant  flux  (the  laser  beam  diffused 
by  a  ground  glass  screen)  for  a  set  of  time  intervals  t^.  The 
negatives  gave  a  set  of  densities  D^,  from  which  the  R-D  curve  was 


obtained  by  setting  log  E  »  log  t  +  k,  with  k  -  log  la  an 

undetermined  constant.  A  cubic  equation  was  found  to  provide  a  good 
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fit  to  this  curve.  Then,  for  other  negatives  on  the  same  film, 
each  D  (hence  each  T)  value  had  a  corresponding  log  E,  obtainable 
to  within  the  unknown  additive  constant  k.  Since  only  relative 

intensities  were  of  interest,  it  was  sufficient  to  set  I  *  KE/t, 
where  the  normalization  factor  K  ■  10  was  a  free  parameter .  The 

procedure  for  selecting  values  of  K  is  described  in  Sec.  2.6.  The 
exposure  time  t  could  be  treated  as  part  of  the  scale  factor  K, 
except  when  background  intensities  were  to  be  subtracted.  In  such 
cases,  t  was  needed  for  photographs  of  the  bubble  scattering  and  of 
the  background  light  to  permit  subtraction  of  a  properly  scaled 
background  intensity.  An  electronic  timer  connected  to  the  shutter 
synchronization  terminal  of  the  camera  measured  t  for  each 
photograph. 

The  transmittance  profile  of  each  negative  was  obtained  as  a 
function  of  absolute  scattering  angle  As  mentioned  previously, 

one  or  more  negatives  on  each  film  contained  a  set  of  exposed  dots 
representing  know  angles  $  .  Microdensitometer  scans  of  these 
negatives  gave  a  series  of  sharp  transmittance  minima  whose  positions 
x^  were  measured.  A  cubic  polynomial  for  $(x)  could  then  be 
obtained  by  least-squares  fitting.  Sprocket-hole  edges  on  the  film 
provided  a  reference  for  the  x  ■  0  position  corresponding  to  the 
edge  of  a  frame.  A  small,  regular  shift  was  detected  in  the 
positions  of  the  sprocket  holes  relative  to  the  camera's  field  of 
view,  accumulating  at  a  rate  of  -9  Um/ frame.  This  was  discovered  by 
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comparing  che  results  of  the  fitted  formula  <t>(x^)  with  the  actual 
angle  data  (^(x^  f°r  films  with  more  than  one  angle  calibration 
frame.  It  became  evident  that  the  difference  ^(x^)  -  cb^Cx^)  had  a 
nearly  linear  dependence  on  the  frame  number  relative  to  the 
beginning  of  the  roll.  This  shift  was  accounted  for  in  the  final 
conversion  of  x  values  into  scattering  angles.  The  camera,  a  Nikon 
F2,  was  mounted  in  such  a  way  as  to  inhibit  angular  motion  which 
might  otherwise  have  occurred  during  the  (motorized)  film 
advancement.  This  angular  shift  appeared  to  be  due  to  the  filling  of 
the  camera's  take-up  reel. 

2 .4 .  Bubble  Shapes 

The  physical-optics  model  of  the  coarse  structure  (Appendix 
A)  assumes  the  bubble  surface  to  be  spherical  at  the  scattering 
plane.  The  degree  of  nonsphericity  of  real  bubbles  is  important, 
therefore,  in  assessing  the  applicability  of  this  model,  and  of  Mie 
theory  as  well.  The  scattering  observed  in  these  experiments  was 
from  near  the  plane  of  the  bubble's  equator;  the  equatorial  plane  of 
a  bubble  is  defined  here  as  the  horizontal  plane  In  which  Its 
diameter  is  a  maximum.  Let  X  denote  the  aspect  ratio  of  a  vertically 
axisymmetric  bubble:  the  ratio  of  its  horizontal  and  vertical 

dimensions.  From  measurements  of  the  equatorial  diameter,  values  of 
X  may  be  estimated  for  the  pendant  and  rising  bubbles  observed. 

For  pendant  bubbled  X ■  a^/z,  where  2a^  is  the  equatorial 
diameter  and  z  is  the  distance  from  the  plane  of  the  equator  to  the 

top  of  the  bubble  (the  bottom  of  the  bubble  is  attached  to  a  needle). 

22 

Numerical  analysis  for  the  shapes  of  stationary  bubbles  allows  z 
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to  be  determined  if  a  is  known.  For  the  largest  pendant  bubble 

P 

observed,  which  had  a^  -  0.981  am,  an  aspect  ratio  of  X  ■  0.969  is 

predicted.  For  smaller  values  of  a  ,  X  approaches  unity.  A  surface 

P 

tension  for  water  of  72.7  dyn/cm  was  assumed  in  using  Qartland  and 
22 

Hartley' s  tables . 

Rising  bubbles  were  trapped  beneath  a  flat  glass  slide  to  be 

measured;  from  the  equatorial  radius  ag  of  the  sessile  (trapped) 

bubble,  the  aspect  ratio  of  the  same  bubble  when  rising  may  be 

estimated.  For  the  largest  sessile  bubble  observed,  which  had  a  * 

3 

22 

0.814  am,  tabulated  values  indicate  a  sphere  of  the  same  volume  has 

radius  a^  *  0.803  am.  Then,  from  the  shape  measurements  made  by 

23 

Siemes  for  air  bubbles  rising  in  water,  an  aspect  ratio  of  X  * 
1.056  may  be  interpolated  for  the  moving  bubble.  Again,  for  smaller 

bubbles  X  tends  toward  unity.  The  equatorial  region  of  observed 

pendant  and  rising  bubbles  should,  by  this  analysis,  present  a  highly 
spherical  scattering  surface. 

The  fine  structure  model  in  Sec.  2.2  assumes  a  bubble  with  a 

circular  cross-section  in  the  scattering  plane;  in  testing  this  model 

the  equatorial  radii  of  bubbles  were  needed.  The  rising  bubble's 

radius  a_  may  differ  somewhat  from  that  of  the  sessile  bubble  a  . 

r  s 

To  estimate  the  difference,  consider  again  the  rising  bubble  with  X  * 

1.056  whose  volume  equals  that  of  a  sphere  with  radius  av  -  0.803 

24 

mm.  Assuming  an  oblate  spheroidal  shape  for  the  rising  bubble 
1/3 

gives  ar  •  ayx  “  0.818  am,  which  differs  from  the  sessile  radius 
ag  ■  0.814  am  by  only  0.5Z.  This  suggests  that  the  error  in 

approximating  ar  by  ag  is  negligible.  The  above  procedure  may 
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slightly  underestimate  the  ratio  ar^as  since  the  water  used  in 

Slemes’  measurements  may  have  been  less  pure  than  ours.  For  rising 

bubbles  X  is  known  to  increase  with  water  purity  for  the  site  range 

23  24 

under  consideration.  * 

Rising  bubbles  were  photographed  «  2  cm  below  the  level  where 

their  size  was  measured.  The  volume  change  due  to  this  height 

difference  may  be  estimated  by  assuming  the  bubble  to  be  an  Ideal  gas 

in  equilibrium  with  the  water.  The  bubble  radius  should  Increase  by 

1/3 

a  factor  (1  +  A P/P)  ,  where  the  change  in  pressure  AP  s  196  Pa 

corresponds  to  a  2-cm  column  of  water,  and  P  -  1  atm  a  10^  Pa;  the 

small  difference  in  heights  makes  this  effect  negligible. 

2.5.  Fine  Structure  in  the  Scattering 

Fine-structure  intensity  oscillations  in  the  scattering  from 
bubbles  may  best  be  observed  near  the  critical  angle  ^ .  Coarse 
oscillations  occur  for  angles  of  <  as  is  evident  in  Fig.  2.1, 

while  as  <t>  continues  to  increase  the  average  intensity  gradually 
decays.  In  the  near  vicinity  of  it)  ,  therefore,  moderately  Intense 

C 

fine- structure  lines  are  visible  without  the  complications  of  a 

superposed  coarse  structure.  For  d)  approaching  180°  the  intensity 

Increases  again  as  "glory"  scattering  phenomena  begin  to  appear;  the 

structures  associated  with  this  region  are  described  in  Chapter  4. 

In  the  present  section,  measured  spacings  of  fine-structure  lines 

near  <j>  will  be  compared  with  the  prediction  of  Eq.  (2.9). 
c 

Figure  2.5  shows  the  experimental  fine-structure  frequency 
(&6)  1  plotted  versus  the  measured  equatorial  radius  for  120  rising 
bubbles  and  23  pendant  bubbles  (note  that  for  rising  bubbles  the 
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PENDANT  BUBBLE  RADIUS  ap  (mm) 


/•s 


RISING  BUBBLE  RADIUS  a,  (mm) 


Figure  2.5.  Measurements  and  model  for  the  angular  frequency 
of  fine-structure  lines.  Data  are  displayed  for  120  rising 
and  23  pendant  bubbles,  whose  radii  as  and  ap,  respectively, 
were  measured  by  microscope.  Mie  results  were  obtained  from 
high-resolution  computations. 
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measured  radius  is  a  ,  as  discussed  in  Sec.  2.4).  The  average 

s 

angular  separation  of  fine-structure  lines  near  $  was  determined 

from  microdensitometer  scans  of  negatives.  The  region  of  measurement 

typically  had  a  width  of  -2.5°. 

The  straight  lines  in  Fig.  2.5  represent  Eq.  (2.9)  adapted  to 

rising  or  pendant  bubbles.  The  values  used  for  and  Xq  are  based 

on  the  average  value  of  the  relative  (air-water)  refractive  index  m 

in  the  experiments.  The  temperature  was  measured  during  the 

25 

experiments,  so  m  could  be  calculated  for  the  laser  wavelength  in 

o 

air  -  6328  A.  Deviations  from  the  average  value  of  a  •  0.75098 

o 

were  negligible,  so  A  *  m  -  4752  A  for  all  the  data.  Also  from 

O  3 

m,  the  critical  angle  <t>c  •  82.649°  is  obtained.  The  corresponding 

incidence  angles  of  the  (0,0)  and  (2,1)  rays  are  0Q  »  48.675°  and  9^ 

■  27.646°,  so  the  sum  of  their  impact  parameters  is  B^  m  1.2150a, 

where  a  is  taken  to  be  a  or  a  ,  for  rising  or  pendant  bubbles, 

s  p 

respectively. 

The  results  displayed  in  Fig.  2.5  compare  favorably  with  the 
fine-structure  model  of  Sec.  2.2.  The  lines  representing  Eq.  (2.9) 
have  a  slope  of  44.625  deg  ^  mm  * .  Linear  least-squares  fits  to  the 
experimental  data,  constrained  to  pass  through  the  origin,  have 
slopes  of  44.444  ±  0.002  deg  *  tun  1  for  the  rising  bubbles  and  45.040 
±  0.035  deg  *  mm  ^  for  the  pendant  bubbles.  In  addition,  six 
fine-structure  measurements  taken  from  Hie  theory  calculations  near 
for  different  bubble  sizes  have  been  included  in  Fig.  2.5;  a 
linear  fit  to  these  results  and  passing  through  the  origin  has  a 
slope  of  44.625  t  0.014  deg  ^  mm  1 ,  in  agreement  with  the 
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aforementioned  result  of  Eq.  (2.9).  The  symbols  used  to  plot  the 
data  in  Fig.  2.5  indicate  the  polarization  of  the  incident  electric 
field  relative  to  the  scattering  plane.  The  effect  of  polarization 
on  the  fine-structure  frequency  is  too  small  to  be  detected  in  these 
experiments,  as  expected  from  Eq.  (2.8)  for  bubbles  of  the  size  range 

a 

considered  here . 

Consideration  of  the  model  leading  to  Eq.  (2.9)  suggests  that 

A<J>  could  be  best  predicted  for  the  rising  bubbles  by  taking  3^  " 

1.2150  ar,  where  a^  is  the  actual  radius  while  rising  as  discussed 

i  in  Sec.  2.4.  Unfortunately  a^  could  not  be  directly  measured; 

however,  the  analysis  of  shapes  in  Sec.  2.4  indicates  that  the 

differences  between  a^  and  ag  should  be  negligible  for  the 

observed  sizes.  The  agreement  with  measured  A<t>  found  by  taking  3^  ■ 

1.2150  a  appears  to  confirm  that  differences  between  a  and  a 
s  r  s 

were  small  for  the  distilled  water  used  in  the  present  experiment . 

Figure  2.6  shows  four  examples  of  scattering  observed  near 

the  critical  angle  from  rising  bubbles  of  different  sizes:  these  are 

positive  reproductions  of  photographs.  The  fine  structure  is  evident 

in  each  case  as  closely  spaced  vertical  lines.  For  Figs.  2.8(a)-(d), 

the  respective  fine-structure  periods  A$  in  degrees  are:  0.1816, 

0.0577,  0.0386,  and  0.4866.  From  Eq.  (2.9)  the  rising  bubble  radius 

af  may  be  calculated.  The  respective  a^  (and  sessile  radii  ag)  in 

mm  are:  0.1235  (0.1171),  0.3890  (0.3880),  0.5808  (0.5817),  and 

0.0461  (0.0436).  The  uncertainties  associated  with  a  measurements 

s 

are  larger  than  those  for  ar  values  derived  from  A<t>  measurements, 


Figure  2.6.  Photographs  of  far-zone  scattering  from  rising  bubbles  in 
water.  Each  photograph  is  aligned  with  the  bottom  scale  showing  the 
scattering  angle  6.  The  bubble  radii  and  incident  polarizations  are 
(a)  ar  -  0.1235  mm,  j  ■  2;  (b)  0.3890  mm,  j  •  2;  (c)  0.5808  mm,  j  -  1; 
(d)  0.0461  an,  j  -  1.  The  corresponding  intensity  profiles  are  shown 
in  Figs.  2.1(b)  and  2. 7-2. 9. 


due  to  difficulties  in  identifying  the  bubble  edges  with  a 
microscope . 

2.6.  Coarse  Structure  in  the  Scattering 

Diffraction  phenomena  are  prominent  in  the  critical-angle 
scattering  region  of  bubbles.  Coarse  oscillations  of  the  scattered 
intensity  are  present  for  followed  by  a  smooth  decay  of 

intensity  as  continues  to  Increase.  Figures  2.6(a)-(c)  each  show 
several  coarse-structure  oscillations,  the  broad,  vertical  light  and 
dark  bands  with  superposed  fine-structure  lines;  they  also  show  the 
gradual  decline  of  average  intensity  at  higher  <fr.  Figure  2.6(d), 
because  of  the  smallness  of  the  bubble,  contains  only  one  faint 
coarse-structure  band,  with  an  intensity  peak  near  75° .  Simple 
geometric  optics,  neglecting  diffraction  effects,  is  unsuccessful  in 
the  critical  region,  predicting  a  cusp  in  the  scattered  Intensity  at 
<t>c,  as  illustrated  in  Fig.  2.1(a).  The  physical-optics  model 
described  in  Appendix  A  is  able  to  approximate  the  coarse  structure 
in  Mie  results  for  a  range  of  bubble  sizes. In  this  section  the 
scattering  observed  from  real  bubbles  will  be  compared  with  this 
physical-optics  model. 

Figures  2.1(b)  and  2.7-2.10  show  experimental  data  for  the 

intensity  as  a  function  of  the  scattering  angle  for  rising  bubbles  of 

different  sizes.  In  each  graph  the  physical-optics  model  result  for 

a  bubble  with  the  same  radius  is  shown  by  a  dashed  curve.  To  compute 

the  model,  the  size  parameter  ka  *  2  ffa/A  is  used,  where  a  Is 

o 

taken  to  be  the  equatorial  radius  ar  of  the  rising  bubble,  a 


quantity  which  was  not  directly  measurable.  Sections  2.4  and  2.5 


SCATTERING  ANGLE  Cd«g> 


Figure  2.7.  Normalized  scattered  intensity  from  a  bubble  with  ka  ■  5144 
and  the  electric  field  parallel  to  the  scattering  plane  (j  •  2  scattering) 
The  solid  curve  is  data  taken  from  the  photograph  in  Fig.  2.6(b),  and  the 
dashed  line  is  the  physical-optics  model. 


SCATTERING  ANGLE  Cdeg3 


Figure  2.8.  Like  Fig.  2.7  but  with  ka  ■  7680,  j  ■  1,  and 
curve  corresponding  to  the  photograph  in  Fig.  2.6(c). 


SCATTERING  ANGLE  Cd«g> 


Figure  2.9.  The  j  ■  1  scattered  intensity  for  ka  »  612.  The  solid 
curve  is  data  from  the  photograph  in  Fig.  2.6(d),  the  dotted  curve  is  the 
Mie  result,  and  the  dashed  curve  is  the  physical-optics  model. 


SCATTERING  ANGLE  <d«g> 


Figure  2.10.  Normalized  scattered  intensity  for  ka  of 
(a)  2748  and  (b)  9699.  These  ka  were  determined  from 
the  respective  radii  of  0.2078  and  0.7335  on  which  were 
inferred  from  the  measured  fine-structure  spacing.  The 
solid  curves  are  measurements;  in  (b)  they  were  smoothed 
to  remove  the  fine  structure.  The  dashed  curve  is  the 
physical-optics  approximation. 


suggest  that  the  sessile  bubble  radius  ag  closely  approximates 
in  these  experiments;  however,  for  computing  the  physical-optics 
model,  *r  was  derived  from  measurements  of  the  fine-structure 
spacing  using  Eq.  (2.9).  Justifications  for  this  procedure  were 

the  following:  (a)  Uncertainties  in  the  derived  a are  smaller 

than  in  the  measured  ag,  due  to  the  difficulty  of  locating  bubble 

edges  by  microscope.  (b)  The  results  of  Sec.  2.5  confirm  that  Eq. 

(2.9)  models  the  relationship  of  Acb  to  the  bubble  radius,  (c)  When 
this  procedure  is  used,  the  model  of  the  coarse  structure  shows 
substantial  agreement  with  the  experimental  data. 

In  the  present  chapter  the  subscript  j  denotes  the 
polarization  of  the  incident  electric  vector  relative  to  the 
scattering  plane;  j  »  1  or  2  for  the  perpendicular  or  parallel 
case,  respectively.  The  model  intensity  is  given  by 


Ij  -  |Sj|2(2/ka)2, 


(2.10) 


where  the  scattering  amplitudes  S j ,  defined  in  Appendix  A,  are 

15  26 

expressed  in  units  commonly  used  '  in  Mie  theory.  The 
normalization  of  the  model  is  chosen  so  that  Ij($)  ■  1  represents 
the  geometric-optics  result  for  a  perfectly  reflecting  sphere®  of 
radius  a.  The  actual  intensity  at  a  distance  R  »  a  from  the 
center  of  the  bubble  is  the  incident  intensity  multiplied  by 
Ij(a/2R)  .  As  discussed  in  Sec.  2.3,  experimental  scattering  results 
taken  from  photographs  were  in  the  form  of  relative  intensity 
profiles  Ir($),  with  a  free  parameter  K  to  be  used  as  a  scale 
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factor  for  Ir.  For  each  of  the  graphs  In  Figs.  2.1(b)  and  2.7-2.10, 

the  constant  K  was  selected  to  optimize  agreement  with  Eq.  (2.10). 

Figures  2.1(b)  and  2. 7-2. 9  contain  the  intensity  profiles 

corresponding  to  the  photographs  in  Figs.  2.6(a)-(d),  respectively; 

from  fine-structure  measurements,  their  respective  bubble  size 

parameters  ka  are:  1633  ±  17,  5144  ±  17,  7680  ±  28,  and  612.0  ± 

2.7.  Additional  intensity  profiles  are  shown  in  Fig.  2.10.  In  each 

case  the  physical-optics  approximation  demonstrates  close  agreement 

with  the  data  in  the  angular  locations  of  coarse  maxima  and  minima. 

Some  differences  between  the  data  and  the  model  are  apparent  in  the 

relative  intensities  of  coarse  oscillations.  These  discrepancies  may 

be  due  In  part  to  the  simplicity  of  the  model,  which  incorporates 

only  two  scattered  rays.  The  model  is  not  an  asymptotic 

2  5  13 

approximation;  as  ka  becomes  large  it  underestimates  ’  ’  the  1^ 

for  4>  >  4>  ,  as  is  evident  in  Figs.  2.7,  2.8,  and  2.10(b).  The 

c 

amplitude  of  coarse-structure  oscillations  is  observed  to  be  greatest 
when  j  ■  2,  as  the  physical-optics  model  pr edicts,  and  the 

fine- structure  amplitudes  are  greatest  when  j  ■  1.  No  figures 

display  both  1^  for  the  same  ka,  but  the  effect  of  j  on  the 
coarse  and  fine  amplitudes  can  be  appreciated  by  comparing  Figs.  2.7 
and  2.8. 


Figures  2.1(b)  and  2.7-2.10  are  representative  examples  of 
the  scattering  data  collected  from  real  bubbles.  Graphs  showing 
similar  agreement  between  the  data  and  the  physical-optics  model  were 
also  obtained  for  bubbles  with  the  following  size  parameters  ka 
(and  radii  ar  in  an):  1364  (0.1031),  1635  (0.1237),  1780  (0.1346), 


2117  (0.1601),  3001  (0.2269),  5899  (0.4461),  6600  (0.4992),  6867 

(0.5193),  7428  (0.5618),  7847  (0.5935),  3995  (0.6803),  and  10  742 

(0.8124).  These  are  presented  as  supplemental  examples  in  Appendix 
C. 

2.7.  Fine-Structure  Contrast  Modulations 

Modulations  in  the  contrast  of  the  fine  structure  were 
usually  visible  in  j  ■  1  scattering  from  bubbles.  These 

modulations  give  some  fine- structure  lines  a  blurred  appearance  In 
photographs;  two  examples  are  conspicuous  in  Fig.  2.6(d)  at  angles  of 
-  74.8  and  77.8°.  Contrast  modulations  were  also  visible  in 
previously  published  Hie  computations:  Fig.  3(a)  of  Ref.  6  and  Fig. 
4  of  Ref.  7.  In  the  present  section  the  angular  period  of  such 
modulations  will  be  approximated,  and  applications  to  bubble  sizing 
will  be  described . 

Consider  a  bubble  with  a  circular  cross-section  of  radius  a 

in  the  scattering  plane;  Fig.  2.2  illustrates  some  of  the  scattered 

rays.  The  far-zone  angular  period  of  Interference  of  two  rays  may  be 

derived  by  the  method  outlined  in  Sec.  2.2.  The  Interference  of  the 

(0,0)  and  (3,1)  rays  has  an  angular  period  of  -  radians,  where 

*  a(sln  8q  +  sin  3^)  is  the  distance  separating  the  two  rays  as 

they-  enter  and  exit  the  bubble.  From  geometric-optics 

considerations,  the  intensity  of  the  (3,1)  ray  taken  alone  is 

expected  to  be  only  about  one- tenth  that  of  the  (2,1)  ray  in  the 

a 

critical-angle  region.  Superposing  the  interference  pattern  of  the 
(0,0)  and  (3,1)  rays  onto  the  fine-structure  pattern  of  the  (0,0)  and 
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(2,1)  rays  should  result  in  contrast  modulations  with  an  approximate 
angular  period  (in  radians)  given  by 

,(A4*)j^  ■  XQ/ (2.11) 
In  the  present  experiments  the  relative  refractive  index  was 

O 

m  ■  0.75098,  so  that  <J)  *  82.649°  and  X  *  4752  A.  The  incidence 

*  c  o 

angles  of  the  (2,1)  and  (3,1)  rays  scattered  to  <}>  are  9„  *  27.646° 

c  2 

and  9^  ■  40.420°,  respectively,  so  B^  -  3^  “  0.1344a.  Comparing 

Eqs.  (2.9)  and  (2.11)  for  this  case  Indicates  there  should  be  -  6.6 
fine-structure  oscillations  within  a  period  (A<t>)^.  The  two  nodes  in 
the  fine-structure  data  of  Fig.  2.9  are  spaced  in  accordance  with 
this  prediction.  In  addition  to  Fig.  2.6(d),  contrast  modulations 
are  visible  in  the  photograph  in  Fig.  2.6(c),  and  in  the  intensity 
profiles  in  Figs.  2.8  and  2.10(a);  these  may  be  used  to  estimate  the 
bubble  size.  There  are  18  modulation  cycles  counted  between  the 
angles  of  80.13  and  84.77  °  in  Fig.  2.8,  giving  ( Ad> ) „  z  0.255°.  From 
Eq.  (2.11)  the  bubble  radius  in  the  scattering  plane  should  be  *• 
0.5790  mm;  this  compares  favorably  with  the  value  of  0.5808  mm 
obtained  from  A4>  in  Sec.  V. 

It  is  possible  to  determine  bubble  sizes  with  high  precision 
by  fitting  observed  intensity  profiles  with  Hie  theory  results.  In 
Fig.  2.9  the  dotted  curve  shows  Mie  results  computed  for  a  bubble 
with  size  parameter  ka  -  612.0;  for  this  case,  the  close  conformity 
to  the  data  is  lost  if  ka  is  changed  in  excess  of  il  in  the  Mie 
calculation.  The  accuracy  of  the  angle  calibration  for  the  data 
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places  an  upper  limit  on  Che  size  for  which  this  fine- structure 

fitting  is  reliable.  The  alignment  of  contrast  modulations  in  Mie 

results  with  those  in  the  data  permits  precise  tea  determinations  up 

to  a  larger  size  limit.  Such  a  procedure  may  become  very  difficult, 

however,  if  ka  is  not  already  known  to  within  ±20.  The  reason  for 

this  is  as  follows.  ?rom  Sq.  (2.3),  the  propagation  phase  difference 

of  the  (2,1)  and  (3,1)  rays  is  ~  02  *  2ka  (3m  cos  p^  -  cos  9^  - 

2m  cos  +  cos  02^’  ^or  scacterii:l8  at  <t>c  with  m  ■  0.75098  this 

becomes  “  ^2  S  ^.161  ka.  If  ka  changes  by  39,  the  phase 

difference  changes  by  2ir,  causing  the  contrast  modulations  at 

<t>c  to  pass  through  one  modulation  cycle.  If  Che  original  uncertainty 

in  ka  exceeds  -  20,  modulation  features  in  Mie  results  may  be 

misaligned  with  those  in  the  data  by  more  than  half  a  period  (A<J>)  . 

M 

2  .8 .  Discussion 

The  physical-optics  approximation  properly  locates  the  angles 

of  coarse  maxima  and  minima.  However,  the  data  obtained  for  large 

bubbles  confirm  that  this  model  underestimates  the  scattered  flux 

into  4>  «  <t>  ,  especially  for  polarization  j  »  1.  Previous 

6  7  12  27 

comparisons  of  '-lie  computations  ’  *  '  with  the  approximation  show 

chat  this  discrepancy  should  depend  not  only  on  ka  and  j,  but  also 

on  the  refractive  index  m.  For  example,  the  mean  intensity  at 

(neglecting  fine  structure)  is  well  approximated  for  bubbles  in 

water*  having  ka  »  100;  but  for  vapor  bubbles  in  liquid  helium  (m  - 

1/1.03,  $c  *  27.7°)  it  is  necessary  to  increase  ka  to  roughly  1000 

07 

to  obtain  similar  agreement.** 
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The  scattering  measurements  are  relevant  to  optical 

techniques  for  detecting ,  sizing ,  and  counting  bubbles  in  vater .  For 

example,  scattering  at  *  90°  is  sometimes  used  for  sizing.^  The 

present  experiments  suggest  that  when  bubbles  are  small  (so  that  It 

is  desirable  to  maximize  the  detected  flux)  it  would  be  preferable  to 

detect  scattering  with  4>  somewhat  less  than  ^  2  32.7°. 

Each  of  the  three  types  of  structure  described  in  Secs. 

2. 5-2. 7  can  be  useful  in  the  sizing  of  bubbles.  The  bubble  radius  a 

in  the  scattering  plane  is  approximately  related  to  the 

fine- structure  period  A<J)  by  Eq.  (2.9)  and  to  the  modulation  period 

(A4>)w  by  Eq.  (2.11).  For  large  bubbles  the  fine  structure  may  be 
M 

poorly  resolved  making  A4>  difficult  to  obtain,  but  ( Ad))^  may  still  be 

measurable.  For  a  wide  range  of  sizes,  use  may  be  made  of 

coarse-structure  oscillations  of  the  model,  which  have  an  angular 

quasi-period  roughly^  ^  -  (  X  /a)^  rad.  Bubbles  may  be  sized  by 

o 

matching  coarse  maxima  and  minima  of  the  data  and  model,  as  in  Figs. 
2.7,  2.8,  and  2.10(b).  Such  a  method  is  possible  even  without  a 
precise  angle  calibration  for  the  data,  since  the  coarse  peaks  near 
<t>c  provide  an  unmistakable  reference  angle.  For  small  sizes  coarse 
oscillations  are  more  difficult  to  use  since  the  range  of  angles 
observed  must  be  large.  For  small  bubbles  it  is  possible  to  obtain 
high-precision  sizing  by  the  alignment  of  fine-structure  peaks  or 
modulations  in  Mle  results  with  corresponding  features  in  the  data, 
as  illustrated  in  Fig.  2.9. 

For  scattering  from  a  polydispersion  of  bubbles  it  is 
expected  that  the  fine  structure  would  be  obscured;  the  coarse 
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structure,  however,  is  less  sensitive  to  bubble  size  and  should  be 
retained  if  the  size  distribution  is  not  too  broad.  Some  photographs 
(not  shown  here)  were  taken  of.  scattering  from  more  than  one  bubble, 
and  these  show  such  effects. 
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APPENDIX  A 


APPROXIMATION  FOR  THE  SCATTERING 


The  main  purpose  of  this  appendix  is  Co  summarize  results  of 

the  physical-optics  approximation^’^  for  the  scattering  amplitudes 

and  intensity  in  the  notation  of  the  present  paper.  Most 

descriptions  of  the  scattering  amplitude  of  spheres  are  phase 

referenced  to  the  point  in  space  corresponding  to  the  center  of  the 

sphere  (see  e.g.  Ref.  15).  That  phase  reference  was  not  used  in 

Refs.  1  and  5;  however,  in  the  present  statement  of  the  phase  [Eqs . 

(A 2)  and  (A6)]  it  will  be  used.  The  final  approximation  for  the 

6  7  12 

intensity  is  identical  to  that  used  in  previous  comparisons  ’  ’  of 
the  model  with  Mie  theory  (for  ka  from  25  to  10  000)  since  the 
intensity  does  not  depend  on  the  phase  reference.  Some  limitations 
of  this  approximation  are  discussed  in  Sec.  2.6. 

The  incident  wave’s  E-field  at  the  point  corresponding  to  the 
center  of  the  sphere  (in  the  absence  of  the  sphere)  is  defined  to  be 
the  real  part  of  5^  exp(-iuit)  •  We  consider  only  cases  where  is 

either  entirely  perpendicular  to  the  scattering  plane  (corresponding 
to  the  assignment  j  *  1)  or  entirely  parallel  to  it  (j  *  2),  since 

these  were  the  polarizations  used  in  the  experiments  and  the  general 
case  may  be  obtained  by  superposition.  The  far-zone  scattered-field 
amplitudes  may  be  written  as  Ej  *  iE^kR)  exp(ikR  -  iujt). 
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Listed  below  are  the  approximate  Sj  which  may  be  inserted  into  Eq. 
(2.10)  to  give  normalized  intensities  Ij  according  to  the  two-term 
physical-optics  model.  Let  S.  ■  S-  .  +  S,  .  where  Sn  .  and  S,  . 
are  amplitudes  associated  with  the  (p,2.)  ■  (0,0)  and  (p,£)  *  (1,0) 
waves,  respectively.  It  is  convenient  to  introduce  where 


spj  '  •ilkaFPj  e,p(lYpj)- 


The  approximation  for  the  reflected  wave  becomes 


■  -2ka  cos  0q  -  H($c  -  <t>)  6 ^ ,  (A2) 

F0j  -  [F(w)  -  F(— )]  2~*  exp(-iir/4),  (A3) 

w  - 

F(w)  a  /  exp(iilT2*)  dz,  (A4) 

0 

w  ■  [(a/X  )  cos  0  ]’ sin(*  -  <!>),  (A5) 

o  c  c 


where  from  Eq.  (2.1),  9^  ■  (  IT  —  <t  )/2.  In  (A2),  5j  is  the  reflection 

phase  shift  given  by  Eq.  (2.5)  and  H  is  a  step  function  giving  H  » 

1  for  <t>  <_<^c  and  H  ■  0  for  <J)  >  An  approximation  to  the 

Fresnel  Integral  F(w)  which  is  useful  for  computation  has  been 
5  13 

previously  cited.  ’  To  appreciate  the  phase  factors  in  (Al)  and 
(A3)  note  that  Fqj  •+■  1  as  w  -*■  ®. 

For  the  single-chord  transmitted  wave  the  approximation^ 


gives 
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■  2ka(mcosp^  -  cos9^), 

FXJ  -  2(1  -  r*)  mc  -  *), 


(A6) 

(A7) 


sin(p1  -  91) 
1  sin(p1  +  9j) 


tan(91  -  p1) 

m  - -  , 

1  tan(91  +  pp 


(A8) 


D  -  i  sin  29,  1 1  -  (m 


-1 


COS  9^008  p^) 


1-1 


/  sin  p. 


(A9) 


For  <j>  >  A  ,  F, 


by 


13 


-  0;  for  Pc»  the  angle  of  incidence  is  given 


tan  9^ 


msin£<t>(l  -  m  cos  $<&) 


-1 


(A10) 


and  the  refraction  angle  by  p^  ■  arcsin(m  3in  0^)  • 

It  should  be  noted  that  the  sign  of  in  Eq*(A8)  differs  from 
that  in  similar  expressions  for  the  Fresnel  coefficients  in 
Eq.(l.lO).  This  is  because  the  particular  sign  convention  used  here 
follows  that  of  the  original  derivation  (Hef.5)  and  was  chosen  to 
give  a  concise  expression  for  the  scattered  fields* 
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APPENDIX  B: 

VIRTUAL  SOURCE  LOCATIONS 

This  appendix  outlines  a  novel  method  for  locating  the 

virtual  source  point  of  a  scattered  ray.  Consider  two  coplanar  (p.J,) 

rays,  R  and  R' ,  leaving  the  spherical  bubble  at  scattering  angles <J> 

and  <J>',  respectively;  the  backward  extrapolations  of  R  and  R' 

will  intersect  at  some  point.  The  virtual  source  point  for  ray 

R  Is  at  this  intersection  in  the  limit  <fr'  *►  To  locate  F  ,  let 

P 

an  x-y  coordinate  system  lie  in  the  scattering  plane  as  shown  in  Fig . 
B.l;the  origin  0  i3  at  the  bubble's  center  and  the  x-axis  is  along 
the  direction  of  the  incident  light.  The  equation  of  the  line 
containing  a  scattered  ray  can  be  expressed  most  simply  by  its  normal 
form.  For  ray  R  this  is 

xcos(<fc  +  tt/2)  +  ysin(<&  +  ir/2)  ■  a(-l)^sin9p,  (Bl) 

where  9p  is  the  angle  of  incidence  (and  departure)  at  the  bubbles's 

surface.  For  R'  an  identical  expression  holds,  except  that  <j>  and 

6  are  replaced  by  4)’  and  9  ' .  These  lines  intersect  at  the  point 
P  P 

F'  whose  coordinates  (x',y')  are  given  by 


(B2) 


x'  »  a(-l)^(sin  0'  cos  $  -  sin  9  cos  )/sin(p  -  p'), 

P  P 

a 

y*  •  a(-l)  sin  9  sec  p  +  x'  tan  $. 

P 

The  coordinates  (Xp,7p)  of  the  virtual  source  point  F^  coincide 
with  (x'  ,y')  in  the  limit  $'-►$,  9p*  -*■  9p.  To  evaluate  (B2)  in 
this  limit,  L 'Hospital's  rule  is  applied  using  the  differential 
operator 


jl  ,  jl  riil  ^£.\  JL 

d6'  "  39'  +  Sq*  +  3p'  d9,;3<J>’  CB3) 

P  P  P  P  P 


3 

39' 


+  2(-ir(pT  -  1) 


3<j>' 


(B4) 


where  (B4)  uses  Eqs.  (2.1)  and  (2.2)  and  the  definition  x  ■  dPp'/dQp' 

»  tan  p  '/tan  9  '.  The  coordinates  of  F  are 
P  P  P 


(d/d9*  )(sin  9'  cos  <t>  -  sin  9'  cos  p') 

lim  _ E _ B - B - 

f 


■4> 

Q'+Q 
P  P 


(d/d9p)(sin($  -  $')] 


y? 


2.  —1 

-a(-l)  sin  9p  sin  p  -  ia(px  -  1)  cos9pCos$, 

+a(-l)^  sin  9p  cos  $  -  ia(px  -  l)"*  cos  9p  sin  6. 


(B5) 


(B6) 

(B7) 


The  exit  plane  of  ray  it  is  normal  to  it  and  tangent  to  the 
bubble;  its  intersection  with  the  scattering  plane  is  shown  in  Fig. 
B«las  the  line  EG,  which  has  the  normal  equation 
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x  cos  <j>  +  y  sin  <J>  ■  a. 


(B8) 


The  poinc  where  R  intersects  its  exit  plane  is  (x_,y  ),  where  from 

D  £ 

(Bl)  and  (B8) 


2, 

xE  -  a  cos  $  -  a(-l)  sin9psin<ti, 

£ 

yE  «  a  sin  4)  +  a(-l)  sin9pcos4>. 


(B9) 

(BIO) 


The  distance  a  from  F  to  the  exit  plane  of  it  is  obtained  by 
P  P 

setting 


ap  "  (xE  "  XF)2  +  (yE  "  yF)2’ 


(Bll) 


which  yields 


ap  -  a[l  +  i(px  -  l)_1cos9p].  (B12) 

This  result  agrees  with  the  previous  analysis  of  the  virtual  source 
location  of  backscattered  glory  waves  (see  Appendix  A  of  Chapter  3). 
It  is  also  applicable  to  the  location  of  the  virtual  line  sources  for 
waves  scattered  by  large  circular  cylinders  illuminated  at  normal 
Incidence  and  to  droplike  objects  as  well. 


SCATTERING  ANGLE  <d«g> 


Figure  C.l.  The  J  •  1  scattered  Intensity  for  lea  ■  13 
The  solid  curve  is  data  froa  a  photograph,  and  the  dashed 
curve  is  the  physical-optics  approximation. 


SCATTERING  ANGLE 


c 


a) 

w 


Figure  C.4.  The  j  ■  2  scattered  intensity  for 
lea  »  2117.  The  solid  curve  is  data  from  a  photograph, 
and  the  dashed  curve  is  the  physical-optics  model. 
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2. 


SCATTERING  ANGLE  <d«s> 


Figure  C.ll.  Like  Fig.  C.l  but  with  ka  *  8995. 
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CHAPTER  3 


FORWARD  SCATTERING 


3.1  Introduction 

The  near-forward  scattered  light  from  a  bubble  contains 
contributions  of  several  types.  This  chapter  is  principally 
concerned  with  the  type  of  scattering  known  as  the  glory;  the  other 
kinds  of  contributions  will  also  be  discussed  briefly.  Photographs 
and  measurements  of  features  related  to  the  forward  glory  of  bubbles 
will  be  presented  for  both  the  near  and  far-zone  scattering .  A 
physical-optics  model  will  be  described  and  its  predictions  will  be 
compared  with  measurements;  a  formally  similar  model  for  the  backward 
glory  was  published  previously  and  is  included  for  reference  in 
Appendix  A.  Some  results  obtained  from  Me  theory  computations  will 
also  be  Included  for  comparisons  with  the  model  and  the  data. 

The  strongest  contribution  to  the  near-forward  scattered 
light  comes  from  what  is  generally  called  forward  diffraction.  For 
any  scattering  object,  that  part  of  the  incident  wavefront  which  is 
not  obstructed  will  produce  a  pattern  of  illumination  extending  into 
the  geometrical  shadow  region  of  the  object.  For  a  large  sphere  of 
radius  a  the  scattering  amplitude  of  the  forward-diffracted  light 
is  approximately  given  by^ 

SppCG)  «  ix2(l  +  cos  0)  J^(x  sin  $)/(x  sin  $), 


(3.1) 
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where  Spp  is  defined  in  Che  same  units  as  Che  amplitudes  of  Eq. 

(l.l).  As  in  Chapter  1,  x  -  2lta/XQ  is  the  size  parameter  of  the 

sphere,  and  $  is  the  scattering  angle  measured  relative  to  the 

forward  direction.  The  3essel  function  J^(z)  is  equal  to  zero  at  z 

»  0,  but  J^(z)/z  »  ^  in  this  limit;  hence  for  exact  forward 

2 

scattering  S_,n(0°)  *  x  /2.  From  Eq.  (1.3)  the  normalized  intensity 
FD 

2 

associated  with  the  diffraction  pattern  i3  proportional  to  x“;  this 
dependence  makes  forward  diffraction  quite  dominant  for  large  values 
of  the  size  parameter  when  4>  is  small.  To  allow  other  types  of 
scattering  to  be  observed  instead,  the  forward-diffracted  light  can 
be  removed  using  polarizers  since  its  polarization  will  be  the  same 
as  that  of  the  incident  light.  This  was  done  in  the  experiments  to 
be  described . 

The  ray  diagram  in  Fig.  3.1  serves  to  illustrate  some  other 

contributions  to  the  near-forward  scattered  light.  The  parameters 

(p,Z)  are  used  to  designate  rays  according  to  their  number  of  chords 

within  the  bubble  and  their  number  of  optic- axis  crossings, 

respectively.  In  Chapter  2  a  physical-optics  model  which  made  use  of 

the  (0,0)  and  (1,0)  rays  was  shown  to  approximate  the  coarse 

intensity  oscillations  observed  near  the  critical-angle  .  The  same 

2  3 

model  approximates  the  coarse  structure  in  Mie  theory  computations  ’ 
for  angles  well  below  However,  as  $  approaches  0s  this  model 

eventually  fails;  new  kinds  of  structure  begin  to  appear  that  are  not 
accounted  for  by  the  interference  of  the  (0,0)  and  (1,0)  rays. 
Though  the  intensities  associated  with  these  rays  do  not  decrease  in 
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the  forward  region,  their  relative  importance  diminishes.  Part  of 

the  reason  for  this  can  be  appreciated  by  considering  the  geometrical 

optics  involved.  For  the  (0,0)  and  (1,0)  rays  the  bubble  acta  like  a 

diverging  lens,  as  Fig.  3.1  illustrates.  The  glory  rays,  on  the 

other  hand,  occurring  both  in  the  backward  and  forward  directions, 

are  scattered  parallel  to  the  optic  axis  but  are  not  coincident  with 

it;  they  are,  in  effect,  weakly  focused  at  infinity.  Because  of  this 

focusing,  it  is  reasonable  to  expect  that  the  intensity  near  the  axis 

at  great  distances  from  the  bubble  would  be  influenced  more  by  the 

glory  scattering  than  by  contibutions  from  diverging  rays.  As 

discussed  in  Sec. 1.2(B),  a  geometrical  divergence  factor  D  i3 

P 

sometimes  useful  in  accounting  for  the  redistribution  of  incident 

energy  by  the  scattering  of  rays.  Finite  values  of  are  found 

for  the  (0,0)  and  (1,0)  cases,  but  for  glory  rays  D  becomes 

P 

infinite  [because  sin  <J»«  0  and  9  *  0  in  Eq.  (1.13)1  •  Geometrical 

P 

optics  cannot  give  the  intensity  at  a  focal  point;  in  the  case  of 

glory  scattering  the  intensity  is  incorrectly  predicted  to  be 

infinite  because  of  the  value  obtained  for  D  .  A  more  sophisticated 

P 

approach,  to  be  detailed  in  the  following  sections,  indicates  that 
the  intensity  associated  with  glory  scattering  is  enhanced  roughly  by 
a  factor  x  (the  size  parameter)  because  of  axial  focusing.  So  the 
near-forward  region  is  expected  to  be  dominated  by  glory  effects  when 
x  is  large  and  the  forward-diffracted  light  is  blocked. 
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3.2.  Physical-Optics  Model  of  the  Forward  Glor\ 


A)  The  Fields  Near  the  Bubble 

It  is  necessary  to  include  interference  and  diffraction 
effects  to  adequately  model  a  -bubble's  far-zone  scattering  close  to 
the  forward  direction.  As  discussed  above,  the  glory  scattering  is 
of  principal  interest  here.  The  physical-optics  approach  to  be  used 
begins  by  representing  the  scattered  light  as  waves  having 
approximately  known  amplitude,  phase,  and  wavefront  curvature  close 
to  the  bubble.  The  present  section  first  examines  these  near-zone 
features  of  the  scattering  which  will  then  be  used  in  formulating  the 
far-zone  model. 

Figure  3.2  illustrates  two  different  glory  rays  and  shows 
some  parameters  important  for  the  wave  description.  An  infinite 
number  of  glory  rays  are  possible,  some  taking  rather  more 
complicated  paths,  but  the  two  shown  turn  out  to  be  the  most 
significant.  The  bubble  is  assumed  to  be  spherical  with  radius  a. 
The  entrance  and  exit  planes  are  tangent  to  the  bubble  and  normal  to 
the  optic  axis.  Along  with  the  (2,0)  glory  ray  in  Fig.  3.2,  two 
other  nearby  (2,0)  rays  are  also  shown.  It  is  evident  from  the 
diagram  that  the  incident  plane  wave  emerges  from  the  bubble  with  a 
curved  wavefront.  The  virtual  source  associated  with  the  (2,0)  glory 


wave  is  at  the  point  F^t  and  the  distance  <x^  is  the  wave’s  radius 
of  curvature  in  the  exit  plane;  a  general  expression  for  for  any 
(p,£)  ray  was  derived  in  Appendix  B  of  Chapter  2.  By  rotating  the 


figure  about  the  optic  axis,  F^  becomes  a  ring  like  source  of  radius 


ENTRANCE 

PLANE 


EXIT 

PLANE 
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b£  and  the  outgoing  glory  wave  becomes  toroidal.  One  of  the 
purposes  of  this  section  is  to  approximate  the  amplitude  distribution 
in  the  exit  plane  due  to  a  glory  wave,  so  that  this  distribution  can 
be  used  later  in  modeling  the  far-zone  scattering  pattern. 
Measurements  will  be  presented  in  Sec. 3. 3  that  support  the 
description  of  the  ringlike  source  for  several  cases. 

The  amplitude  of  a  glory  wave  leaving  the  bubble  depends  on 
losses  due  to  transmissions  and  internal  reflections,  and  also  on 
changes  in  the  area  over  which  the  incident  energy  is  eventually 
spread.  Let  the  incident  light  be  linearly  polarized  with  amplitude 

Ei  exp(-iu)t).  It  is  useful  to  define  a  pair  of  orthogonal  basis 

*  A 
vectors  e^  (h  *  1,2)  in  the  entrance  and  exit  planes:  let  e^  be 

parallel  to  the  direction  of  the  incident  wave's  polarization  and  % 

be  perpendicular  to  it.  The  amplitude  of  the  glory  wave  in  the  exit 

plane  will  be  given  by  its  components  with  respect  to  these  vectors. 

Points  in  the  exit  plane  are  specified  by  the  polar  coordinates 

(s.tIO,  where  s  is  the  radial  distance  from  C'  (where  the  optic  axis 

intersects  the  exit  plane),  and  $  is  the  pola»  angle  measured  from 

the  direction  of  e^.  The  amplitude  of  the  (p,2.)  glory  wave  in  the 

i  A  2 » 

exit  plane  is  given  by  E  ■  E  e,  +  E  e,  where 

P  Pi  P  2 

Ep  "  Eiqp*  ^  exPtinp  +  bp)2/2ap5»  (3.2) 


where  q  characterizes 
P 

includes  the  transmission 


the  spreading  of  the  wavefront  and  F*1 
and  reflection  effects.  The  exponential 
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term  gives  the  approximate  phase  at  points  in  the  exit  plane,  as  will 
be  discussed  below.  An  expression  similar  in  form  to  Eq.  (3.2) 
appears  in  Eq.  (1)  of  Appendix  A  where  the  baclcward  glory  is 
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desrlbed;  the  factor  q^  ■  a^/(  ap  ~  *s  c^e  8ame  the  present 

ll 

case.  The  amplitude  factors  F  are  derived  using  van  de  Hulst’s 

method  of  decomposing  the  incident  light  into  perpendicular  and 

4 

parallel  components  with  respect  to  the  scattering  plane  .  The 
resulting  forms  are 


12  2 
F  »  sin  cos 

F2  ■  i(c2  -  Cj)  sin  2xb» 


(3.3) 


where  the  coefficients  c^  and  c^  characterize  the  effects  of 
transmission  and  reflection  on  the  perpendicular  (j  ■  1)  and  parallel 
(j  ■  2)  components: 


c  . 
J 


(i  -  r2). 


(3.4) 


The  F*1  then  express  the  resulting  amplitudes  with  respect  to  the 
more  convenient  basis  vectors  8,^.  The  r^  in  Eq.  (3.4)  are 
Fresnel's  coefficients  as  given  in  Eq.  (1.10),  evaluated  using  Che 
incidence  and  refraction  angles  9p  and  for  the  glory  ray  of 
interest.  In  general,  these  angles  are  determined  numerically  from 
Sqs.(1.8)  and  (1.9),  but  exact  trigonometric  solutions  are  known  for 
several  cases,  including  thirteen  forward  glory  rays;  these  solutions 
are  presented  in  Appendix  B  of  this  chapter.  The  sign  factor  in  Eq. 
(3.4)  differs  from  the  corresponding  expression  in  Appendix  A  because 
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Che  forward  glory  does  not  exhibit  a  geometrical  inversion  of 

amplitudes  present  in  the  backward  case.  An  interesting  angular 

2 

dependence  can  be  seen  from  Eq.  (3.3)  for  F  ,  the  amplitude  factor 
for  the  light  having  a  polarization  perpendicular  to  that  of  the 
incident  light.  This  "cross-polarized"  scattering  is  predicted  to 
have  zero  amplitude  at  angles  of  \J)  *  0  ®  ,  ±90°  ,  and  180®  in  the  exit 
plane,  with  measured  from  the  direction  of  the  incident  light's 
polarization.  Photographs  will  be  presented  in  the  following 
sections  of  this  chapter  that  illustrate  a  similar  angular  dependence 
at  the  ringlike  sources  and  also  in  the  far-zone  scattering. 

The  exponential  expression  in  Eq.  (3.2)  approximates  the 
phase  of  the  scattered  field  at  the  radial  distance  s  in  the  exit 
plane.  The  first  term, 

np  -  2ka(l  -  cos  9p  +  mp  cos  pp)  -  $u(p  +  2  -  1),  '  (3.5) 

accounts  for  the  total  phase  delay  of  the  (p,  2)  glory  ray  as  it 

propagates  from  the  entrance  plane  to  the  exit  plane,  including  phase 

changes  due  to  focal-line  crossings;  these  phase  considerations  were 

2 

discussed  in  Sec. 2. 2.  The  other  term,  k(s  -  b  )  /2  a  ,  gives  an 

P  P 

additional  phase  delay  for  s  *  bp.  It  is  obtained  by  approximating 
the  toroidal  wavefront  as  a  quadratic  surface  with  radius  of 
curvature  at  the  exit  plane. 

The  amplitude  and  phase  factors  described  above  are  to  be 
evaluated  for  each  glory-wave  of  Interest.  The  total  field  In  the 
exit  plane  due  to  glory  scattering  is  then  approximated  by  summing 


the  Ep  from  Eq.  (3.2).  But  the  f ar-zone  approximation  which  follows 
is  most  conveniently  carried  out  by  using  the  individual  glory-wave 
amplitudes  and  then  summing  the  resulting  fields. 

3)  The  F ar-Zone  Scattering 

The  glory-scattered  E-field  at  a  distant  point  Q  close  to 
the  forward  axis  will  be  approximated  using  Che  amplitude 
distribution  in  the  exit  plane,  Eq.  (3.2).  Let  Q  be  at  a  distance 
R  from  the  point  C',  where  the  exit  plane  is  tangent  to  the  bubble 
on  the  optic  axis.  Figure  3.3  shows  the  coordinate  system  that  will 
be  used.  The  scattering  angle  4  to  point  Q  is  assumed  to  be  small, 
while  R  is  very  large;  these  parameters  are  shown  out  of  proportion 
in  the  figure. 

The  field  at  Q  is  obtained  using  far-zone  diffraction 
theory  by  integrating  the  contributions  at  Q  from  all  points  (s.'l)  ) 
in  the  exit  plane,  which  is  show  as  the  plane  x'-y’ .  The  procedures 
for  performing  the  integration  have  been  given  in  several  papers; 
Appendix  A.  outlines  the  method,  and  detailed  discussions  are  provided 
in  Refs.  3  and  6.  At  Q  the  field  due  to  the  (p , i)  glory  wave  is 

^  i  2  A 

approximately  E^  ■  *1  +  Ep  *2  wh,re 

E*  -  (k/2iriR)  E*  ql*  D*  exp(ikR  +  in  ),  (3.6) 

P  1  P  P  P 


with 


-  W*1  bp(Aoap)^exp[iir/4-  ikap(l  -  co»$)] 


(3.7) 
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and 

W1  •  ir[(c,  +  c-)Jn(u)  +  (c.  -  c5)  J?(u)  cos  2C] 

1  2  0  1  2  2  (3.8) 

VT  -  ir(Cl  -  c2)  J2(u)  sin  25. 

The  argument  of  the  Bessel  functions  Jq  and  J 2  is  u  *  kb^  sin  <|>  s 
and  5  is  the  azimuthal  angle  of  point  Q,  as  shown  in  Fig.  3.3.  Like 
the  exit-plane  polar  angle  i|> ,  the  far-zone  azimuthal  angle  5  is 
measured  with  respect  to  the  direction  of  the  incident  polarization, 
which  was  defined  by  e^.  And  like  the  near-zone  amplitudes,  the 
E-fleld  at  Q  is  expressed  by  its  components  parallel  (h  *  1)  and 

perpendicular  (h  »  2)  to  the  direction  of  the  incident  polarization. 

2 

From  examination  of  Eq.  (3.8)  for  the  amplitude  factor  W  , 

2 

it  is  evident  that  the  far-zone  field  Ep  vanishes  at  angles  5  * 
0°,  ±90°,  and  130°.  A  dark  cross  is  thus  predicted  in  the  h  ■  2 

scattering,  with  arms  parallel  and  perpendicular  to  the  direction  of 
the  incident  polarization.  The  terminology  used  here  refers  to  the 
h  «  2  scattering  as  "cross-polarized" ,  not  because  of  this 

distinctive  scattering  pattern  but  because  this  light  has  a 

polarization  which  is  transverse  to  that  of  the  Incident  light,  in 
the  sense  described  in  Ref .7.  In  the  experiments  to  be  described  it 
was  the  h  ■  2  glory  scattering  that  was  observed  in  the  far  zone, 
because  the  forward-diffracted  light  was  removed  by  a  technique  using 
crossed  polarizers.  Equation  (3.8)  also  predicts  a  ring  pattern 
centered  on  the  forward  axis  for  both  polarizations.  In  the  h  -  2 

case  the  rings  are  expected  to  be  spaced  according  to  the  period  of 
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Che  J2  Bessel  funtion,  while  for  Che  h  •  1  case  oscillaCions  of 
boch  Che  Jq  and  funccions  are  involved. 

The  far-zone  model  ouclined  above  may  Include  as  many  (p,Z  ) 
glory  waves  as  are  desired.  The  cocal  approximace  field  ac  Q  due 


Co  glory  waves  is  obcained  by  summing  Che  Ep  from  Eq.  (3.6)  for 
all  che  (p,£)  cases  included: 


gh 

glory 


(R,4>,§) 


l 

(P.A)  P 


(3.9) 


The  relaCive  incensicy  of  Che  glory  scactering  ac  Q  i3  Chen  simply 


I 


h 

glory 


"glory 1 


(3.10) 


IC  is  inscruccive  to  compare  the  intensities  which  result 
when  each  glory  wave  is  considered  alone.  From  Eqs.(3.6)  -  (3.8), 
the  maximum  value  of  the  h  ■  1  intensity  from  a  single  (p,£)  glory 
wave  is 


IpC<0  -  0a)  -  2xIR(bp/a)2(^^)7r(Cl  +  c2)2, 


(3.11) 


since  Che  Jq  Bessel  funecion  has  its  maximum  value  of  1  at  u  *  0. 
The  maximum  h  ■  2  intensity  arising  from  a  single  (p,£)  glory  wave 
is  given  by 


a -a 


lj(*  -  *g)  -  2xIR  (bp/a)2(^2-I)7r[0.4865(c1  -  c2)]2, 


(3.12) 
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where  the  azimuthal  angle  5  has  been  set  to  45°  and  $  ■  arcsin 

s 

(3.0542/kb  )  Is  the  scattering  angle  at  which  the  first  maximum  of 
P 

O 

the  Bessel  function  occurs:  J2(3.0542)  z  0.4865.  In  Eqs. 

(3.11)  and  (3.12)  the  normalization  factor  1^  Is  the  Intensity 

scattered  by  a  perfectly  reflecting  sphere  of  radius  a,  as  given  by 

h 

Eq.  (1.2).  Table  3.1  compares  the  normalized  I  for  a  number  of 

glory  waves  each  taken  separately.  The  size  parameter  x  *  4000  was 

used  with  each  of  two  relative  refractive  index  values.  For  m  * 

(1.403)  ^  this  corresponds  to  a  bubble  of  radius  a  s  287  um  in  a 

silicone  liquid  that  was  used  In  experiments  to  be  described.  The 

same  value  of  x  but  with  m  «  0.75  corresponds  to  a  bubble  of 

radius  a  a  302 y m  in  water  (taking  the  light  to  have  a  wavelength  in 

air  of  632.8  nm  for  both  examples).  Some  of  the  focal  parameters 

associated  with  each  of  the  glory  waves  are  also  Included  in  the 

table.  It  should  be  emphasized  that  the  intensities  I  h  are  not  to 

P 

be  summed  as  a  means  of  obtaining  the  total  glory-scattering 

intensity.  The  amplitudes  must  be  summed  as  in  Eq.  (3.9).  3ut  Table 

3.1  serves  to  illustrate  some  qualities  of  the  individual  terms.  It 

demonstrates  that  some  of  the  (p,£)  glory  waves  when  taken  separately 

can  yield  an  intensity  that  exceeds  the  Intensity  from  a  perfectly 

reflecting  sphere  by  a  considerable  factor.  It  also  suggests  that 

the  main  contributions  to  the  total  glory  scattering  must  come  from 

those  waves  having  small  values  of  the  parameters  p  and  £  . 

Figure  3.4  shows  the  physical-optics  model  for  the  far-zone 
2 

intensity  I  in  comparison  with  the  results  of  the  exact  Mie 
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Table  3.1.  Model  Results  for  x  ■  4000  with  a  ■  (1.403)“**  and  a  ■  0.75 


■ 

C P'l ) 

•» 

(deg) 

b/a 

a/a 

* 

(6  -  0*) 

$ 

(4  •  4g) 

(1.403)"1 

(2,0) 

40.30 

0.647 

1.093 

224.8 

11.97 

(3,0) 

43.37 

0.687 

1.035 

41.21 

1.749 

(4,0) 

44.32 

0.699 

1.019 

12.56 

0.507 

(5,0) 

44.74 

0.704 

1.012 

5.048 

0.200 

(6,0) 

44.96 

0.707 

1.008 

^ 

0.094 

(7,0) 

45.10 

0.708 

1.006 

1.292 

0.050 

(4,2) 

19.06 

0.327 

1.095 

0.006 

0.006 

(5,2) 

28.28 

0.474 

1.061 

0.004 

0.001 

0.75 

(2.0) 

43.57 

0.689 

1.093 

278.5 

9.571 

(3,0) 

46.57 

0.726 

1.034 

48.36 

1.402 

(4,0) 

47.49 

0.737 

1.018 

14.54 

0.405 

(5,0) 

47.89 

0.742 

1.011 

5.813 

0.159 

(6,0). 

48.11 

0.744 

1.008 

2.766 

0.075 

(7.0) 

48.24 

0.746 

1.006 

1.481 

0.040 

(4,2) 

20.32 

0.347 

1.101 

0.003 

0.003 
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Figure  3.4.  Normalized  cross-polarized  intensities  as  a  function  of 
scattering  angle.  The  dotted  curve  is  a  physical-optics  model  using 
only  the  (2,0)  glory  wave;  the  dashed  curve  is  the  model  result 
incorporating  the  (2,0)  and  (3,0)  waves;  the  solid  curve  is  the  Mie 
theory  result. 


theory.  The  size  parameter  is  x  ■  4000,  and  the  refractive  index  is 
m  *  0.75.  The  case,  which  is  dominated  by  forward  diffraction, 

is  not  included.  The  dotted  curve  in  this  figure  is  the  model  result 
when  only  the  (2,0)  glory  wave  is  used;  its  maximum  intensity  agrees 
with  the  value  of  9.571  listed  in  Table  3.1.  The  dashed  curve  is  the 
physical-optics  model  incorporating  both  the  (2,0)  and  (3,0) 
contributions.  The  solid  curve  is  the  Mie  result,  normalized  to  I 
*  1  like  the  other  two  curves;  this  normalization  of  Mie  results  was 
discussed  in  Sec.  1.2(A).  The  Mie  theoretic  expression  for  the 

cross-polarized  intensity  is  given  by  Eqs.  (6)  and  (7)  of  Ref.  7. 
Evidently,  the  inclusion  of  the  (3,0)  wave  gives  an  improved 
approximation  for  the  intensity  as  compared  with  the  Mie  cheory.  The 
angular  locations  of  the  maxima  and  minima  are  well-approximated  by 
the  simpler  model  using  only  the  (2,0)  wave;  the  incorporation  of  the 
(3,0)  term  gives  no  appreciable  improvement  here. 

9 

The  Mie  theory  curve  in  Fig.  3.4  is  similar  to  one  published 
for  a  spherical  bubble  in  water  with  x  »  3040.  These  and  other  Mie 
computations  confirm  that  for  spherical  bubbles  in  water  the 
cross-polarized  near-forward  scattering  should  be  quasi- per iodic  in 
for  a  wide  range  of  bubble  sizes. 

3.3  Observations  of  Virtual  Ringlike  Sources 


The  experimental  arrangement  for  observing  forward  glory 
effects  from  bubbles  is  diagrammed  in  Fig.  3.5.  The  basic  design  was 
the  same  for  the  near-  and  far-zone  observations,  except  that 
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different  camera  lens  systems  were  employed  in  the  two  cases.  The 
incident  light  was  from  a  15  mW  He-Ne  laser.  This  beam  was  passed 
through  a  polarizing  prism  to  make  it  highly  plane-polarized .  A 
glass  cell  containing  the  host  liquid  for  the  bubble  was  immediately 
followed  by  a  second  polarizing  prism.  This  allowed  the  polarization 
of  the  observed  scattered  light  to  be  selected.  Photographs  of  the 
scattering  were  made  by  a  camera  focused  at  the  desired  range.  This 
section  describes  measurements  of  the  glory-wave  virtual  sources  that 
appear  within  the  bubble.  Photographs  of  these  focal  circles  were 
taken  using  a  Nikon  bellows  extender  and  a  reversed  50  mm  Nikkor  lens 
to  allow  short-range  focusing  with  about  5X  magnification. 

The  bubble  was  injected  with  a  syringe  into  a  silicone  oil 
(Dow  Corning  200  Fluid)  of  high  kinematic  viscosity  (  -  600  000 
centistokes) ,  so  that  it  remained  almost  immobilized  but  assumed  a 
nearly  spherical  shape.  The  liquid  had  a  refractive  index  of  1.403. 
3ubbles  were  produced  with  diameters  ranging  from  -  1  to  7  mm.  The 
upper  limit  corresponded  to  the  maximum  laser  beam  diameter  allowed 
by  the  apertures  of  the  polarizing  prisms. 

Photographs  of  the  virtual  ringlike  sources  inside  a  bubble 
are  shown  in  Fig.  3.6.  The  bubble  radius  was  a  ■  1.375  mm.  Parts 
(a)  and  (b)  show,  respectively,  the  co-polarized  (h  *  1)  and 
cross-polarized  (h  ■  2)  appearances  of  the  rings.  For  both 
photographs  the  direction  of  the  incident  polarization  was  vertical. 
Recall  that  the  physical-optics  model  for  the  h  ■  2  exit-plane 
amplitude  predicted  zeros  at  polar  angles  of  0°,±90°,  and  130°  [Sq. 
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(a) 


Figure  3.6. 
Photographs  of  the 

(a)  co-polarized  and 

(b)  cross-polarized 
forward  glory  circles 
for  a  bubble  with  a 
radius  of  1.88  nan  in 

a  liquid  of  refractive 
index  1.403. 


(b) 
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(3.3)] •  The  ring  sources  in  Fig.  3.6(b)  manifest  four  dark  regions 
with  this  same  kind  of  angular  dependence. 

It  is  possible  to  distinguish  several  different-sized  rings 
in  the  negatives  from  which  Fig.  3.6  is  taken.  Microscope 
measurements  of  these  focal-circle  radii  are  listed  in  Table  3.2  with 
similar  measurements  from  other  bubbles.  From  geometrical 
considerations,  the  virtual  ring  source  of  a  (p,2)  glory  wave  is 
expected  to  have  a  radius  b^  ■  a  sin  0^,  where  8^  is  the  (p,2)  glory 

ray's  angle  of  incidence.  The  values  listed  in  the  table  give  the 

ratio  of  the  ring  radius  to  the  bubble  radius  for  the  visible  rings 
counted  from  smallest  to  largest.  The  same  data  is  presented  in 

graphical  form  in  Fig.  3.7.  The  first  (smallest)  ring  can  be 
associated  with  the  (2,0)  ray.  It  is  considerably  more  Intense  than 
the  others;  because  of  this  it  appears  broader  in  the  photographs  and 
its  exact  radius  was  more  difficult  to  determine.  More  precise 
measurements  were  able  to  be  made  on  the  second  ring;  its  radius  fit 
the  theoretical  b/a  prediction  for  the  (3,0)  glory  wave  source. 
The  higher  order  rings  became  quite  dim  while  their  spacing 

decreased,  making  measurements  of  their  radii  less  precise  again. 

The  horizontal  lines  in  Fig.  3.7  are  the  theoretical  values  of  b/a. 

P 

Their  spacing  illustrates  the  way  the  focal  circles  coalesce  as  the 
value  of  p  increases. 
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Tabla  3.2.  Focal-Circle  Radii  for  Bubbles  in  Silicon*  Oil 

a  bj/a  b3/a  b4/a  bg/a 

(am)  (1st  Ring)  (2nd  Ring)  (3rd  Ring)  (4th  Ring) 


0.533 

1.290 

1.823 

1.875 

1.993 

2.046 

3.501 


0.645  ±  0.006 
0.651  t  0.004 
0.652  ±  0.004 
0.651  ±  0.004 
0.649  ±  0.003 
0.655  ±  0.010 
0.656  ±  0.006 


0.684  ±  0.004 
0.688  ±  0.004 
0.686  ±  0.004 
0.683  ±  0.004 
0.689  ±  0.009 
0.688  ±  0.002 


0.696  ±  0.003 
0.699  ±  0.002 
0.694  ±  0.004 
0.696  ±  0.003 
0.700  ±  0.008 
0.701  ±  0.002 


0.707  ±  0.004 
0.701  ±  0.004 
0.707  ±  0.004 


Theory  0.6468 


0.6867 


0.6986 


0.7039 


Ring  Radius  /  a 
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3.4  Observations  of  Far-Zone  Scattering 

k)  Single  Bubbles  in  Silicone  Oil 

The  far-zone  scattering  from  single  bubbles  was  observed 
using  the  experimental  arrangement  shown  In  Fig.  3.5.  The  camera  had 
a  200  tan- focal- length  lens,  which  was  focused  at  *  to  allow  the 
far-zone  scattering  pattern  to  be  recorded.  The  bubbles  were  In 
silicone  oil  having  refractive  index  1.403.  Their  diameters  ranged 
from  -  0.2  to  1.9  mm,  while  the  laser  beam  diameter  could  be  adjusted 
to  as  large  as  7  mm.  The  beam  was  made  to  be  highly  collimated  in 
the  vicinity  of  the  bubbles  by  using  a  parallel-plate  shearing 
interferometer.10  The  polarizing  prisms  used  were  of  elllpsometric 
grade;  attempts  to  observe  the  far-zone  forward  glory  with  polarizers 
of  lesser  quality  were  unsuccessful  because  of  a  large  transmitted 
background.  The  scattering  angles  that  could  be  observed  were 
limited  to  ♦  $4°  by  the  aperture  of  the  polarizer. 

Figures  3.8  and  3.9  show  the  cross-polarized  (h  ■  2)  far-zone 
patterns  for  bubbles  of  radius  a  *  0.203  tan  and  a  *  0.600mm, 
respectively.  The  most  distinctive  features  of  the  photographs  are 
the  four  dark  lobes  that  form  a  cross  pattern  like  that  predicted  by 
Eq.  (3.8).  The  direction  of  the  incident  polarization  was  vertical 
for  both  photographs. 

The  angular  spacing  of  the  dark  rings  was  measured  for 
these  and  several  other  photographs,  to  compare  them  with  the  ring 
structure  predicted  In  Eq.  (3.8).  For  small  angles  $,  u  s  kb  4>  ; 
after  the  first  few  oscillations  the  zeros  of  the  Jj  Bessel  function 
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1  dej 


Figure  3.8.  Photograph  of  the  cross-polarized  near-forward  scattering 
from  a  bubble  of  radius  0.203  ran  in  silicone  oil. 


Figure  3.9.  Photograph  of  the  cross-polarized  near-forward  scattering 
from  a  bubble  of  radius  0.600  oa  in  silicone  oil. 
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are  spaced  by  approximately  it,  so  A<t>  z  ir/kb.  Taking  b  ■  b.„:  m 
0.6468a  gives 

A$  s  1.998  x  10"2(deg/mm)  a"1  (3.13) 

as  the  model  prediction  for  the  approximate  angular  spacing  of  the 
rings.  Figure  3. 10 shows  the  average  measured  ring  spacing s  from  11 
photographs  like  Figs.  3.8  and  3.9.  Also  included  on  the  graph  are 
measurements  taken  from  Mie  theory  calculations  at  3  bubble  sizes. 
The  measurements  show  substantial  agreement  with  Eq.  (3.13),  though 
this  prediction  is  a  rather  simplistic  one  since  it  takes  only  the 
(2,0)  glory  wave  into  account.  It  is  important  to  note  that  the  ring 
structures  in  the  far-zone  glory  scattering  have  a  considerably  wider 
spacing  than  would  rings  associated  with  the  forward  diffraction 
pattern.  From  Eq.  (3.1)  it  is  clear  that  the  diffraction  rings  would 
have  an  angular  spacing 


(AdOpp  s  ir/(ka),  (3.14) 

which  is  smaller  by  a  factor  b/a  than  the  glory-ring  spacing  found 
in  Eq.  (3.13). 

Figure  3.9  shows  evidence  of  interference  effects  between 
waves  in  the  near-forward  direction;  the  ring  structure  is  modulated 
in  intensity  at  certain  scattering  angles.  It  has  not  been 


determined  whether  these  modulations  arise  from  interference  of  the 
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•  Photographs 
■  Mie  Theory 


Figure  3.10.  Angular  spacing  of  dark  rings  in  the  far-zone 
cross-polarized  scattering.  Measurements  were  made  from 
photographs  like  Figs.  3.8  and  3.9  and  also  from  Mie  theory 
computations.  The  line  is  the  physical-optics  model 
prediction  in  Eq.  (3.13). 


separate  glory  waves  or  from  the  presence  of  (0,0)  and  (1,0) 
scattered  waves. 
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3)  Clouds  of  Bubbles  in  Water 

The  forward-glory  scattering  from  clouds  of  small  bubbles  in 
water  was  photographed  using  the  apparatus  of  Fig.  3.5  with  slight 
modification.  The  glass  cell  was  fitted  with  a  tungsten  wire 
stretched  across  near  the  bottom,  and  this  was  used  in  forming 
bubbles  by  electrolysis.  The  other  electrode  was  a  needle  Inserted 
into  the  water  at  the  top  of  Che  cell.  The  camera  had  a  200 
mm- focal- length  lens  focused  at  infinity. 

Figure  3.11  is  an  example  of  Che  cross-polarized  forward 
glory  scattering  observed.  The  photograph  shows  the  distinctive  dark 
cross  pattern  that  is  predicted  by  the  physical-optics  model  for  each 
bubble  individually.  The  first  few  dark  rings  are  also  visible.  It 
is  likely  that  a  dispersion  of  bubble  sizes  was  present  and  tended  to 
obscure  the  ring  structure  beyond  these  first  intensity  minima.  No 
direct  measurements  were  made  to  cry  to  determine  an  average  bubble 
3ize,  but  an  estimate  may  be  made  from  the  spacing  of  the  observable 
rings.  On  the  negative  from  which  Fig.  3.11  was  obtained,  the  ring 
spacing  is  about  1  mm.  With  a  focal  distance  of  200  on  this  gives  an 
angular  separation  of  about  0.005  rad.  From  Eq.  (3.13)  applied  to 
the  case  of  water  (where  b^  •  0.6893)  the  radius  a  is  estimated  at 
about  69  urn,  which  is  a  reasonable  result.  From  this  observation  it 
can  be  concluded  chat  freely  rising  bubbles  of  this  approximate 
radius  and  smaller  are  sufficiently  spherical  for  the  forward  glory 
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Figure  3.11.  Cross-polarized  near-forward  scattering  from  a 
cloud  of  small  bubbles  rising  in  water. 


to  out.  Very  largo  bubbles,  on  thd  other  band,  take  on 
.oo.bh.rid  .hope.  ..  diseased  to  Sec. 2.4.  It  i.  to  be  ..petted 
that  the  symmetry  of  the  cro.s-polarized  near-forward  scattering  fro. 
large  bobbles  will  differ  from  that  seen  in  fig.  3.11. 
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Obeervtdoas  of  llfht  beokscsttared  from  air  bubbler  la  a  viaeoua  liquid  dauuxuirau 
an  enhancement  dua  co  axial  focusing.  A  phyilcsl-opaea  appro xlmartoa  for  tha  eroaa- 
polarlxad  aaaaarlaf  correctly  daacrlbaa  tha  tpastag  of  regular  faaruraa  observed.  Tha 
aos-anaa  polarixad  seactarlag  la  wx  adequately  daaortbad  by  a  ilagla  elaaa  of  rays. 


PACS  aumbara:  12.20.0g,  12. 10. Ho,  92.10.Pt 


Tha  Mia  solution1  for  electromagnetic  scatter- 
tng  by  a  sphere  frequently  doas  not  land  to  diraet 
interpretation  of  tha  angular  scattering  partam. 
Consequently,  modal*  hara  baan  developed  to  fa- 
eilitata  an  undarataoding  at  tha  structura  In  tha 
scattarad  intensity  prasaot  where  tntansuy  ia 
ploctad  aa  a  function  of  tha  seattarlng  angla  o  or 
tha  sixa  parameter  x  «*a  (A  Is  tha  wars  nuabar; 
a  1*  tha  sphara  radius).  Thaaa  modal*  hara  am- 
phaslsad  tha  angular  ragiona  whara  diffraction  1* 
important  for  a  drop  of  water  in  air:  tha  rain* 
to»,w  <r»  NOV  and  tr*  0*.»  »  In  tha  seattarlng 


of  light  by  s  spherical  sir  bubble  in  a  liquid  or  in 
glass,  tha  raal  part  of  tha  rafractir*  indax  of  the 
spherq  la  lass  than  that  of  the  surroundings  and 
tha  models  must  ba  significantly  modified.  New 
phenomena  appear,  such  as  diffraction7''  in  the 
region  of  tha  critical  scattering  angle  <dc.  Hera 
we  report  the  first  detailed  observation*  of  back- 
seattartng  by  air  bubbles  In  liquid*  and  give  a 
modal  which  describes  some  of  the  observed  fea¬ 
tures.  We  refer  to  this  aa  glory  because,  a*  In 
the  ease  of  drops, ^  tha  <i »  130°  seattarlng  is  en¬ 
hanced  when  x  1*  large. 
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Van  da  Bulat1'4  pn  a  partial  explanation  a i  tba 
enhancement  (or  drops  by  noting  tba  axial  focusing 
of  thoae  bacfcacsttored  raya  which  hava  a  nonzero 
impact  parameter.  When  modeling  thla  focusing 
in  the  tar  (laid,  diffraction  provides  an  essential 
correction  to  ray  optica  because  the  (actor  in  the 
scattered  intensity  which  accounts  (or  geometri¬ 
cal  divergence  of  the  raya  goes  to  *  ae  uj  -  180' . 
Examination  of  this  (actor  in  ray-optics  models 
at  scattering  by  bubbles'  shows  that  thla  infinity 
is  not  restricted  to  drops.  We  have  modeled  the 
bac  It  scattering  with  a  physical-optics  approxima¬ 
tion.  The  procedure  la  to  (a)  compute  amplitudes 
in  an  exit  plane  In  contact  with  the  bubble  via  ray 
optica,  and  lb)  allow  this  wave  to  diffract  to  the 
(ar  (laid  where  the  distance  (ram  the  bubble’s 
center  S  »>  Ha1. 

Figure  1  illustrates  several  raya  which  lead  to 
backaeatterlng.  The  paths  are  determined  by  the 
number  of  chorda  o  and  m  "a,/*,,  where  the  re¬ 
tractive  indices  at  the  Inner  and  outer  media,  m , 
and  hi  ,,  ire  taken  to  be  real.  Figure  1  is  drawn 
with  m ' 1  >1.403  which  corresponds  to  an  air  bub¬ 
ble  in  the  dimethyl- tUaxane- polymer  liquid  used 
in  the  experiment.  All  raya  satisfy  sin#  attxi. 
For  <c  ■  180" ,  tba  aff-axls  (or  glory )  raya  hava  a 
•9  and  pod,  whare40  *2-p)90*,  g  la  a 

nonnagxtlve  integer  (g  « o  (or  raya  in  Fig.  1),  and 

*  t  requires  that »  *  J.  The  exit  plane  (dashed 
Une  in  Fig.  V  touches  C'  with  its  normal  parallal 
to  the  propagation  direction  at  the  incident  wave. 

Our  description  at  the  field  In  the  exit  plane  is 
facilitated  by  considering  the  propagation  of  a 


no.  l.  Itaya  which  aeatrlburt  to  baskacattarmr. 
The  local  angle  o(  incidence  is  *  and  C  is  the  bubble’s 

cantor. 


wavelet  da  which  Ilea  cloaa  to  tha  hackee  altered 
path.  Figure  1  shows  dt  for  8  *3;  it  emerges  as 
curve  Ve‘ .  This  curve  appears  to  coma  from  a 
ringlike  source  at  F  known  aa  tha  focal  circla  ia 
tha  analogous  8*2  scattering  from  drops*  with 

<  m  <  2.  Tha  source  is  ringlike  because  the 
figure  may  be  rotated  around  the  CC'  axis.  The 
radius  of  the  ring  la  t>  *a  sind.  After  the  incident 
ray  crosses  the  dashed  vertical  plane  (the  en¬ 
trance  plane),  tha  propagation  phase  delay  for 
reaching  tha  exit  plane  is  1  -hall  -con*  -(1 
- cosd laectt  -J)  »2in8 eoao  1.  The  ray  croasea 
the  exit  plane  at  a  radius  a  from  C'  with  s/a  -  Had 
-(1  -co*J)tan<8  -J).  The  radlua  a  of  arc  d'f 
fallows  from  the  curvature  at  j  -S:  a  **4fV 
da*)'1  *«{t  f^(pr  -  O' 1  eoai],  where  r  -tan}/ 
tank.  The  spreading  of  ttwwavelst  la  eharactar- 
lzad  by  g  *Umfcf’«"ie)  as 3*-0  where  the  bar  de¬ 
note*  the  are  length.  An  equivalent  expression 
(or  ?  islllmlh  -*<0)1/ (8  -a  aim))!  asd-J;  its 
value  from  L 'Hospital's  rale  ia  a/ (a  -a).  Vectors 
i,  (1 » 1,2)  denote  orthogonal  basts  vectors  in  both 
tha  entrance  end  the  salt  planes;  r,  is  chosen 
parallal  to  the  polarization  at  the  incident  wave's 
electric  (laid  £,  «xp<- tu<). 

In  tha  exit  plus,  the  (laid  B,‘i ,  of  the  outgoing 
8th  glory  ware  is  computed  by  applying  Van  de 
Hulst’s  method  ed  first  decomposing  the  fields 
perpendicular  and  parallel  to  the  scattering 
plane.1'4  Exit-plane  polar  coordinate*  centered 
an  C’  are  (*,»),  where  *  ia  the  angle  relative  to 
»i  and  i  and  i  denote  local  basis  vectors.  Wa  as¬ 
ms*  i*l  and  in  Franal't  coafflctems  r,  (or 
tha  internal  rtfleettoM  whare  i  *1,2  for  (laid* 
parallal  to  *  and  I,  respectively.  tit-4i«a, 
the  multiple  internal  reflections  give 

£,'  *£fg" ''*F‘  exptiD  *I*(S  i,  (1) 

where  8«t*  *  n<jn  »/D,  F‘(»)  «c,  auA  -c.coe1*, 

£*(*)  •*<«, -c,) sin2-,  and  c,  *(- 'M{1 
The  new  phase  term  *  accounts  for  tha 
crossing  of  caustics  or  'focal  fine*";  its  value 
la1'**  -»(p*g)/2.  Tha  r,  are  evaluatad  at  b. 
r^rtati *4)  tad  r,  small  -5)/taa<J 
*i).  Tha  sign  factor  la  c,  accounts  (or  a  geo¬ 
metrical  inversion  (preseat  whan /  *2  and  p  is 
odd)  wMch  is  not  evident  in  descriptions  of  p  s2 
glory  la  drop*.**4 

The  field  at  a  distant  pah*  9  ti  computed  aa 
follow*.  Tha  left  extension  at  the  CC'  axle  makes 
aa  angle  v  with  C’if.  Whm>  la  email  and 
•£'»*#*,  scalar  diffraction  theory  and  the  Fraun- 
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holer  approximation1*  give 

e  ) 

e'  i&tk-p7* 

*  (,*sW'axp(i*(s -d)'/2a  Ids,  (2) 

'•</'>(,  F‘«xpi-<*«  star  coa<»  -oWit,  (3) 

where  4  is  the  ingle  between  #,  and  the  projection 
of  CQ  oo  the  exit  plane.  In  Eq.  (2),  the  approxl- 
mation  given  by  Eq.  (i)  has  been  extended  beyond 
its  useful  domain  in  anticipation  of  the  stationery 
phase  approximation  (SPA)  at  the  integral.  Street 
evaluation  of  Eq.  (3)  gives  >W»(y,4,s  *4) 

•e{(e,  *ciV,<n)  *(ct  -  c,y,fc)  coa24l  and  **  •r(el 
-c,y,(e)  stn24,  where  u  >**  liny.  The  SPA  al 
Eq.  (2)  gives  the  nth  glory  contribution  to  the  scat¬ 
tered  field  when  *0J  o ,  and  thus  x ,  are  large,  la 
the  experiments  to  be  described  x  a  4000  and  the 
SPA  is  applicable. 

The  total  (laid  may  be  approximated  by  summing 
the  £,'  (ram  Eq.  (2)  with  the  (is Ids  due  to  axial 
reflections  and  surface  waves.  Surface  wave  con¬ 
tributions  should  be  small  (or  the  observed  bub¬ 
bles  because  of  the  largeness  of  x.  To  determine 
which  glory  and  axial  terms  are  Important  to  the 
total  (laid,  and  (or  other  beurlacic  reasons,  con¬ 
sider  the  (-polarized  intensity  of  the  pth  Held 
tahm  alont.  The  SPA  at  Eq.  (2)  gives 

f,,e(2/r)x/i/;i,[lt’,(y,4jj*,  (4) 

where  /,  ef.aV 4Ea  la  the  total  intensity  at  a  dis¬ 
tance  A  •£§  from  a  prrfrclty  rtfltcting  sphere 
of  radius  a  predicted  by  ray  optics.*  I,  is  the  In¬ 
cident  lotaoslty,  and/,,<**1a/aJg>41(a-a)/«». 

In  Eq.  (4),  A  has  replaced  A'  (ram  (2)  and  v  be¬ 
comes  180“  -vi  because  A  »o.  Geometrical  op¬ 
tics'''  gives  the  Intensities  /,'  of  separate  axial 


reflections  (a.g.,  t>  *0  and  2  In  fig.  1)  which  are 
proportional  to  a*.  The  strongest  reflection  has 
p  >0  and  /  ■  l;  (or  y  >0,  */,4»i  -  1)*/&h  ♦  i)* 

while  /,'  -0.  Sines  doss  not  dapend  on  a,  /,' 

■*  *aJ  and  glory  terms  dominate  the  becfcacattering 
when  u  la  large. 

Consider  a  bubble  with  r  •*000  and  <»  >1.403*  *. 
The  atrongset  glory  terms  have  g  >0  and  p  >3,  4, 
and  S;  the  /,'.  I,  lory  >0  are,  respectively,  1.03, 
0.43,  and  0.16.  The  /,*  decrease  with  increasing 
p  as  a  result  of  the  partial  reflections  in  the  bub¬ 
ble.  The  strongest  axial  ray  gives  >0.028. 
The  Interference  of  the  fields  depends  on  a  and 
our  Mia  computations  verify  that  the  oackacat- 
tered  Intensity  in  not  simply  proportional  to  a* 
even  (or  this  large  value  of  x .  The  (  >2  (cross-- 
polarized)  scattering  is,  however,  nearly  domi¬ 
nated  by  then  >3  glory  term.  Beeanse  of  sym¬ 
metry,  /  >2  scattering  vanishes  as  r-0.  The 
/,'(y  >0,4)  have  maxima  at  4  ■*  43*  and  *  133°  and 
they  vanish  at  4  >0",  *  90°,  and  180*.  Lety-y, 
locate  the  first  maxima  of /,'(y,‘ v4J°).  The 
largest  I  >3  terms  hsve  4  *43°)//,  >0.33  and 
OJO  (or  P  >3  and  4.  To  the  extent  that  p  >3  scat¬ 
tering  may  be  neglected,  the  (  >2  intensity  will  be 
quamperiodie  lay. 

We  have  numerically  verified  the  validity  of 
Eq.  (4)  by  uateg  Debye’s  localisation  principle1'* 
to  modify  Mie  theory  so  chat  only  partial  waves 
associated  with  t>  >3  rays  were  included  in  the 
Mis  series.  Furthermore,  when  Eq.  (4)  is  ap¬ 
plied  to  spheres  with  certain  m  >  1,  the  resulting 
/,'(r  >4)  agree  wtth  the  glory  ’analog  '  tabulated 
is  Ref.  11.  This  analog  was  derived  by  applying 
the  Watson  treoWormatton  to  tha  y  >0  Mis  series. 

rigor e  2  diagrams  the  experiment.  A  syringe 
injected  bubbles  into  the  liquid.  The  liquid  had  a 
high  kinematic  viscosity  i»  600000  cS;  l  stabs  (S) 
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Voiumi  47.  Sum 
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no.  4.  MonraMcn  aad  modal  lor  Ilia  angular 
saparstkm  of  tba  dark  flags  la  tba  1  *  2  acattarlac. 


-1  cm1/ sac  |  aad  a  single  bubbla  could  ba  obaarrad 
(or  hours.  Tha  laaar’s  pcarar  output  was  5  aW 
and  tba  baam  dtamatar  was  3am.  Tha  warwlanglh 
in  tba  liquid,  23 /k,  was  (832.8  na)/l.t03;  #,  lay 
In  tba  spiittar’s  plana  at  Ineidooca.  Tha  eamarm 
was  tocuasd  on  -  so  that  tba  photographs  racord- 
ad  tha  far-fleld  intanaity  pattern.'- **  Photocrap ba 
war*  aada  with  a  »  0.3 -0.3  aa  corraapondlaf  to 
x  *  4000-11000.  Expoaur*  Haas  war*  typicaUy 
S  s  (or  TrUC  (11m  and  a  200-mm-(ocal-L*nftb 
caaara  Iona. 

Fleur*  3  damonatrates  that  tba  acattarlac  baa 
roufhlr  tba  dapandanc*  on  4  pradktad  by  Eq.  (4); 

4  •O'1  corresponds  to  seatterlnc  toward  tba  top  at 
tha  photographs  and  y  «0*  corraaponda  to  tba  can¬ 
tor  at  tba  synmacry.  Fleur*  Mb!  shows  that  tha 
l  - 1  seatterlnc  (or  v  *  0 X  la  sigmflatmly  strong¬ 
er  (or  4  ••  90”  than  it  la  (or  (  *0*.  Thin  as raoa 
with  tha  foUmrtng  modal  r* suits:  (1)  (c^c,)**  l 
((or  a  *3  w*  predict  cyc,»  -3  .2);  and  (11)  (or  this 
x ,  tha  /,*  dapood  only  weakly  an  4  and  are  domi¬ 
nated  by  the  /,*,  On*  prediction  at  Eq.  (4)  could 
b*  quant! tartrely  cbackad:  whan  both  aim'  ■  y  and 
a  » 1,  tba  minima  in  should  ba  spaced  by  Ay 
rad  such  that  Wa»"i,  where  (or  *  «3,  i/a 
•0.447.  Fleur*  4  comparts  this  with  the  mean 


T—4~t  at  >40  dark  rlncs  lylnc  outald*  tha  9th 
rlnc  (ram  tha  canter.  The  error  bar*  combine 
uncertainties  in  measured  a  aad  Ay  with  tboaa  at 
earxwettans  dn*  to  ratractlaa  at  tha  call-air  Inter¬ 
lace'  aad  tba  tilt  at  tba  call.  Fleur*  4  shows  that 
9  >3  rays  ■*"■»<—«»  tba  t  *2  acattarlac.  Tha  modu¬ 
lation*  at  tba  Intensity  alone  4  **  43*  la  Fig.  3(b) 
show  that  other  ray*  contribute  to  ( •  1  acattarlac 
sine*  tba  predicted  *  [,/0M  1*. 

lb  conclusion,  backacatterlnc  tram  bubbles  can 
ba  enhanced  by  axial  focusing.  Tha  number  at 
sleadlcaac  glory  terms  dapanda  on  m .  The  mala 
eomributlans  differ  (ram  those  (or  water  drops 
where  surface  wares*  and  other  diffraction  re¬ 
lated  terms'  play  aa  aaaantlal  rote,  a  focusing 
wars  not  prasant,  seatterlnc  by  large  bubbles 
would  bo  «  fa  la  tho  region7-*  (<pe  -10*) » <g  i  180*, 
whar*  <3,  *2  coo"  ‘n  “  89*  (or  <*■'  -1.403.  W*  also 
(lad  arl dance  at  p  «3  glory  In  Mia  computations 
for  bubbla*  In  .rotar. 
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APPENDIX  B 


INCIDENCE  ANGLES  OF  SCATTERED  RAYS 

This  appendix  provides  some  solutions  for  the  incidence 
angles  of  scattered  (p,£)  rays  when  the  scattering  angle  <t>  and  the 

relative  refractive  index  ra  are  known.  In  general,  9^  can  be  found 

numerically  from  Eqs.  (1*9)  and  (1.9),  but  a  precise  determination 
may  require  numerous  iterations.  A  few  exact  trigonometric  solutions 
for  9p(<t>,m)  have  been  found,  and  these  are  given  in  the  following 
tables.  For  arbitrary  values  of  $ ,  9p  may  be  computed  for  (p,2)  rays 
having  p  £  3  using  the  procedures  given  in  Table  B.l.  For  <t>  »  0  °  , 

solutions  of  0p  for  thirteen  different  (p,Z)  cases  are  given  in  Table 

B.2.  In  both  tables,  the  solution  gives  the  value  of  sin  9 

P 

Ellison  and  Peetz  (Ref.  11,  p.118)  give  expressions  for  9  for  the 

P 

forward  (2,0)  and  (3,0)  rays.  While  their  expression  for  the  (2,0) 
ray  is  equivalent  to  the  one  given  here,  their  result  for  the  (3,0) 
ray  is  incorrect,  apparently  due  to  a  misprint. 


Table  B.l.  Incidence  Angles  of  (p.l)  gars 


(p  ,1) 

sin  0_  ■  b  /  a 

P  P 

(0,0) 

cos  & 

(1.0) 

2  .1 
o  sin  ( 1  +  nr  -  2  m  cos  i<fc)”7 

(2,0)  , 
(2.1)  1 

(16  -  R2)^R/8,  where: 

R  -  A  +  Q  -  [8  +  3A2  -  2m2  -  Q2  -  2AQ~1(4  ♦  m2  -  A2)] 
Q  -  [A2  +  2(4  -  m2)(l  +  cos  Y)/3]^, 

Y  ■  3“ 1  arccos  [27  m4  (1  +■  cos$)(4-  m2)-3  -  1], 

A  »  2“^  m[cos£$  +  (-l)^sinH], 

(3,0) 

(3.1)  } 

(3.2) 

3R  -  4R3,  where: 

R  -  (3  -  A  -  QcosY)*,  Q  -  (4A2  -  6A  -  3  m2  +  9)*, 

Y  -  3""1  arccos  {Q-3[(3 -A)3  -  9(3  +  A)(m2  -  A2)]}, 
A  -  aeMK-l/ct/e-  ff/12)  +  (3  -  2£)ir/12]. 
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Table  B.2.  Incidence  Angles  of  (p,&)  Forward  Glorr  Rara 

(p,A)  sin  8p  -  bp  /  a 
(2,0)  in2[l  +  (1  +  8  a"2)*] 

(3,0)  *n3/2(l  +  3n-1)* 

[tf2]  >  «[Q  *  (i  -  Q2  +  AQ-1  «/32)i'l.  * 

where:  Q  -  6“*[1  -  (2.5)*  cos  Y]*, 

T  -  3-1  arccos  [-(2.5)“*(1.4  +  27n2/320)]. 

(5!°)  }  i®[2.5  +  A(1.25  +  An)*]*  * 

(7,0)  in3~*(7  +  7*  cos  Y  +  21*  sin  Y)* 

(7.2)  in3"*(7  +  7*  cosY  -  21*  sin  Y)* 

(7.4)  in3”*(7  -  28*  cos  Y)*, 

where  T  ■  3”*  arccos  [i  7-3^2(7  -  27Ao)]. 

(9,0)  in  [9  -  2Q  +  (27  -  4Q2  +  3Q"1)*]* 

(9.2)  in  [9  +  2Q  +  (27  -  4Q2  -  3Q"1)*]* 

(9.4)  in  (9  +  2Q  -  (27  -  4Q2  -  3Q“X)^]^ 

(9,6)  in  [9  -  2Q  -  (27  -  4Q2  +  3<f  *)*]*, 

where:  Q  -  i[9  -  8(3  -  4A  n/3)*  cos  Y]*.  # 

Y  -  3"1  arccos  [-4.5(1  -  An) (3  -  4An/3)"3/2 

*  In  «ch  of  tl»a.  con..,  *  -  (*\  |"  £  I  I  or  ll 
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CHAPTER  4 


DISCUSSION 


The  light  scattering  properties  of  bubbles  have  been  examined 
theoretically  and  experimentally  with  emphasis  on  three  scattering 
regions:  the  critical  angle,  and  the  forward  and  backward  directions. 
In  each  of  these  angular  regions  the  scattering  pattern  is  not 
descrlbable  by  the  methods  of  simple  geometrical  optics  because  of 
the  prominence  of  diffraction  and  interference  phenomena.  If  the 
bubble  can  be  assumed  to  be  a  homogeneous  dielectric  sphere,  an  exact 
solution  to  the  scattering  problem  is  provided  by  the  Ilie  theory. 
Its  results  are  valuable,  but  the  Mie  theory  has  certain  drawbacks  as 
well.  It  is  computationally  intricate  and  time-consuming,  and  it 
gives  no  direct  Insight  into  the  physical  origin  of  features  in  the 
scattering .  Physical-optics  approximations  have  been  developed  which 
can  predict  and  explain  noteworthy  features  in  the  scattering ,  and 
are  more  readily  adaptable  to  changes  in  the  size,  shape,  and 
refractive  index  of  the  bubble.  When  dealing  with  less  idealized 
scattering  situations  the  physical  insight  may  be  of  considerable 
value. 

The  scattering  features  which  have  been  photographed  and  modeled 
are  among  the  most  outstanding  to  be  observed  from  bubbles;  the 
scattering  in  each  of  these  diffraction  regions  is  notably  bright. 
The  transition  to  total  reflection  of  light  near  the  critical  angle 
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Is  a  consequence  of  Che  relative  refractive  Index  of  the  bubble  being 
less  than  that  of  the  surrounding  liquid.  The  enhanced  intensity  of 
the  forward  and  backward  glory  arises  due  to  the  spherical  symmetry 
of  the  bubble.  The  physical-optics  approximations  given  in  this 
dissertation  allow  the  scattering  in  these  transition  regions  to  be 
modeled,  and  provide  useful  techniques  for  the  sizing  and  detection 
of  bubbles.  Scattering  features  such  as  the  fine  structure  near  the 
critical  angle  can  be  used  for  high-precision  sizing  of  single 
bubbles,  while  Che  coarse  structure  in  that  region  gives  information 
about  average  sizes  in  a  distribution  of  bubbles.  The  ring 
structures  in  the  forward  and  backward  glory  scattering  are  a  measure 
of  bubble  size  and  are  also  sensitive  to  effects  of  nonsphericity. 
The  dark  cross  pattern  associated  with  the  glory  scattering  provides 
a  distinctive  means  of  detecting  the  presence  of  bubbles,  whether 
singly  or  in  an  ensemble. 


